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NOTATION 


Vectors  and  scalars  are  generally  represented  by  lower  case 

letters  and  matrices  by  upper  case  letters.  The  derivative  with  respect 

to  time  of  a  vector  x  is  denoted  by  x.  The  transpose  of  a  matrix  A  is 
T 

denoted  by  A'*'.  The  absolute  value  of  a  scalar  f  is  denoted  |f|.  Norm  of 
a  matrix  or  vector  M  la  denoted  BMII  .  Notation  for  a  set  V  composed  of  all 
objects  u  having  some  property  P  is 

V  =  {u  :  u  has  property  P}. 

k  jL 

A  function  f(p)  is  0(vi  )  if  there  exist  positive  constants  y  and  c  such  that 
the  norm  fl  fll  satisfies 


ifn  <  cy^ 

when 

k>  0 

1  f  1  >  cy^ 

when 

k<  0 

for  all  ye  (0,y*].  The  expression  Q>0  means  that  Q  is  a  positive  definite 
matrix.  The  i-th  eigenvalue  of  a  matrix  R  is  denoted  as  X  (R) . 


1.  INTRODUCTION 


!.!♦  Problem  Description 

In  the  last  twenty-five  years  adaptive  systems  theory  has  received 
a  lot  of  attention  by  theoreticians  and  practitioners.  Global  stability  of 
MRAS  (Model  Reference  Adaptive  Schemes),  an  open  problem  for  almost  two 
de’cades,  was  recently  solved  for  both  continuous  and  discrete  SISO  (Single- 
Input  Single-Output)  systems  [1-16].  Moreover,  different  adaptive  schemes 
developed  Independently  have  been  shown  to  be  equivalent  and  diverse  approaches 
have  been  unified  [13],  [17-19].  However,  there  still  remains  a  significant 
gap  between  the  available  theoretical  methodologies  and  potential  applications 
of  such  adaptive  schemes.  Global  stability  properties  are  guaranteed  under 
the  "matching  assumption"  that  the  model  order  Is  not  lower  than  the  order  of 
the  unknown  plant.  Since  this  restrictive  assumption  Is  likely  to  be  violated 
In  applications.  It  Is  Important  to  determine  the  robustness  and  performance 
of  adaptive  schemes  with  respect  to  such  modeling  errors. 

Recently  several  attempts  have  been  made  to  formulate  and  analyze 
reduced-order  MRAS  [20-24].  The  results  of  such  studies  depend  on  th^  ,  - — 
characterization  of  the  model-plant  mismatch. 'l!hlqf  thesis  Is  divided  Into  six 
sections.  Different  characterizations  are  given  for  the  model-plant  mismatch 

ru 

In  Sections  2  to  5  and  6.  -60^ characterization  In  Sections  2  to  5  assumes  a 
separation  of  time  scales  between  the  modeled  and  unmodeled  phenomena.  The 
order  of  the  model  Is  equal  to  the  order  of  the  slow  part  of  the  unknown 

plant  and  the  model-plant  mismatch  Is  due  to  the  fast  part  of  the  plant.  In 

p 

most  applications  the  slow  part  consists  of  "dominant^  modes,  while  the 
neglected  fast  modes  are  considered  as  parasltlcs.  Ovf  approach  Is 


asymptotic.  \  A  crucial  parameter  In  our  analysis  Is  the  "speed  ratio"  u  of 
the  slow  versus  the  fast  phenomena.  The  limit  as  p->0  means  that  the  fast 
part  of  the  plant  reaches  Its  steady-state  Instantaneously,  that  Is,  the 
plant  reduces  to  Its  slow  part.  The  fictitious  "reduced-order"  plant  Is 
thus  obtained  when  In  the  actual  plant  y  >  0  Is  replaced  by  y>0.  This 
singular  perturbation  approach  Is  a  convenient  parameterization  of  the  model- 
plant  mismatch  because  It  allows  an  asymptotic  analysis  using  the  limit  as 
y-^0.  In  our  formulation  adaptive  schemes  are  designed  for  the  reduced  order 
plants  and  then  applied  to  the  actual  plants.  They  are  considered  robust  If 
the  error  In  their  performance,  due  to  model-plant  mismatch.  Is  0(y).  In 
Section  6  the  model-plant  mismatch  Is  characterized  by  neglected  Interconnec¬ 
tions  between  different  subsystems.  Adaptive  schemes  are  first  designed  for 
the  decoupled  subsystems  and  then  applied  to  the  overall  system  with  the 
Interconnections.  ^ 

A 

1.2.  Section  Review 

In  Section  2  we  analyze  the  behavior  of  continuous- time  Identifiers 
and  adaptive  observers  with  respect  to  modeling  errors,  consisting  of  fast 
parasltlcs  which  are  weakly  observable  In  the  plant  output.  The  adaptive 
schemes  considered  are  shown  to  be  robust  provided  the  Input  signal  Is 
dominantly  rich,  that  Is  It  Is  rich  for  the  dominant  modes,  but  does  not 
contain  high  frequencies  In  the  parasitic  range.  The  bounds  Indicate  possi¬ 
bilities  for  reducing  the  error  by  a  proper  choice  of  the  Input  signal. 

The  performance  of  reduced-order  discrete-time  parallel  and  series 
parallel  adaptive  Identifiers  with  respect  to  weakly  observable  parasltlcs  is 
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treated  in  Section  3.  While  in  a  deterministic  environment  with  no  modeling 
error  the  two  schemes  give  Identical  results,  their  performance  is  different 
when  used  in  a  deterministic  environment  with  modeling  error.  Detailed 
bounds  on  the  composite  output/parameter  error  are  established.  The 
dependence  of  these  bounds  on  the  initial  identification  error,  the  charac¬ 
teristics  of  the  input,  and  the  speed  ratio  u  are  shown  to  be  crucial. 

In  Section  4,  it  is  shown  that  the  weak  observability  assumption 
used  in  Sections  2  and  3  is  crucial.  That  is  when  the  parasltlcs  are  strongly 
observable  the  adaptive  schemes  are  no  longer  robust.  A  redesign  procedure 
has  been  Introduced  which  re-establishes  the  robustness  of  identifiers  and 
adaptive  observers  with  respect  to  fast  strongly  observable  parasltlcs. 

The  effects  of  unmodeled  fast  dynamics  on  the  stability  and  per¬ 
formance  of  adaptive  control  schemes  are  analyzed  in  Section  5.  In  the 
regulation  problem  global  stability  properties  are  no  longer  guaranteed,  but 
a  region  of  attraction  exists  for  exact  adaptive  regulation.  In  the  case  of 
tracking  a  more  robust  adaptive  law  is  proposed  which  guarantees  the  existence 
of  a  region  of  attraction  from  which  all  signals  converge  to  a  residual  set 
which  contains  the  equlllbrl\im  for  exact  tracking.  The  size  of  this  set 
depends  on  design  parameters,  the  frequency  range  of  parasltlcs,  and  the 
reference  input  signal  characteristics. 

In  Section  6  the  problem  of  regulation  and  tracking  of  a  class  of 
large-scale  linear  time  invariant  systems  with  unknown  parameters  is 
considered.  An  approach  is  developed  for  adaptive  regulation  and  tracking 
using  decentralized  adaptive  controllers.  Sufficient  conditions  are  estab¬ 
lished  in  a  form  of  algebraic  criteria  which  can  guarantee  stability  or 
boundedness  under  certain  structural  perturbations. 
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2.  CONTINUOUS-TIME  IDENTIFIERS  AND  ADAPTIVE  OBSERVERS 

2.1.  Characterization  of  the  Mismatch 

Suppose  that  an  (n-fm)th  order  linear  time  invariant  plant  has  n 
slow  (''dominant")  and  m  fast  ("parasitic")  modes,  that  is  n  of  its  eigenvalues 
are  0(1)  and  the  remaining  m  are  0(-^),  where  u  is  a  small  positive  scalar. 
Without  loss  of  generality  such  a  plant  can  be  represented  by  the  standard 
singular  perturbation  form 

i  »  A^j^x  +  A^2*£  (2.1) 

WXj  ■  A^J^X  +  +  B^u  (2.2) 

where  x,x^  are  n  and  m  vectors,  respectively,  and  u  is  an  r  input  vector. 

In  this  form  the  dominant  and  the  fast  parts  do  not  appear  explicitly,  but 
state  variables  x,  x^  and  parameter  u  have  clear  physical  meaning.  Typically 
y  represents  small  time  constants,  masses,  inertias  etc.  [25]  .  It  is  known 
from  [25]  that  the  contribution  of  the  fast  modes  to  x  is  only  0(u)  and  hence 
X  can  be  used  as  the  state  of  the  dominant  part  of  the  plant.  State  x^  is 
formed  of  a  "fast  transient"  and  a  "quasi-steady  state"  defined  as  the 
solution  of  (2.2)  with  yx^*0.  This  motivates  the  definition  of  the  fast 
parasitic  state  as 

n  *  +  Lx  +  A^H^u  (2.3) 

xriiere  L  is  required  to  satisfy  [26 ] 

A22L-A2^  +  yLAj^2^ -  yLAj^j^  -  0.  (2.4) 

From  (2.4)  and  assuming  that  A22  exists  we  see  that 

L  -  A22A2j^  +  0(y) 


(2.5) 


and  hence,  for  y  small,  n  Is  the  difference  between  and  Its  "quasi-steady 
state"  -A22A2j^x- A22B2U.  Defining 

A  ■  A^^”A^2L,  B  ■  "*  ^12^f  ^f  *  ^f  "  ^22  ^'^^*^12* 

-  B2  +  wLBj^,  H  -  Aj^2  (2.6; 

and  substituting  (2.3)  Into  (2.1),  (2.2)  we  obtain  a  representation  of  (2.1) 
(2.2)  with  the  dominant  and  parasitic  parts  appearing  explicitly 


Ax  +  Bu  +  Hn 

(2.7 

A^n  +  yA^H^u 

(2.8 

as  shown  In  Fig.  2.1.  The  parasitic  state  n  enters  the  dominant  part  as  a 
disturbance  Input  Hn.  The  parasitic  part  of  the  plant  Is  driven  by  the 
Input  uu.  The  steady-state  gain  In  the  disturbance  path  Is 

-HA“Sj  (2.9] 

and  will  be  shown  to  play  a  crucial  role  In  the  Identification/observation 
error  due  to  parasltics.  From  (2.7),  (2.8)  and  Fig.  2.1  It  can  also  be 
anticipated  that  another  crucial  factor  will  be  the  fequency  content  of  u(t) . 
If  u(t)  contains  frequencies  as  high  as  0(-^),  then  yu  Is  0(1)  and  the  steady 
state  of  n  Is  0(1),  that  Is  the  Influence  of  the  parasltics  Is  significant. 

We  can  now  formulate  the  model-plant  mismatch  problem  as  follows . 

I 

The  goal  Is  to  Identify  and  observe  only  the  nth  order  dominant  part  of  the 
plant.  Hence  the  order  of  the  adjustable  model  is  chosen  to  be  n.  An 
adaptive  scheme  Is  then  designed  Ignoring  the  parasltics,  that  is  assuming 
y«0.  Our  problem  Is  to  analyze  the  performance  of  this  scheme  when  applied 
to  the  plant  with  parasltics,  that  Is  with  y>0. 


We  address  this  problem  under  the  following  assumptions : 
Assimptlon  I;  The  plant  with  parasltlcs  Is  asymptotically  stable,  that  Is 


ReX(A)<0,  ReX(A^)<0. 


(2.1( 


Assumption  II;  The  pair  (A,B)  of  the  dominant  part  Is  controllable. 

Noting  from  (2.5), (2. 6)  that  A,  A^,  and  B  are  differentiable 
functions  of  u,  for  u  sufficiently  small  assumptions  I  and  II  will  hold  If 
A^  Is  replaced  by  and  A  and  B  are  replaced  by 


“  ^12^22® 


-1, 


'll  “  *12*22^21 


'l  “12"22"2 


(2.i: 


which  are  the  matrices  of  the  plant  without  parasltlcs , 


X  -  AqX  +  BqU  (2.12 

obtained  by  setting  y  -  0  In  (2.7),  (2.8)  and  In  Fig.  2.1. 

Assumption  III;  Input  u(t)  and  Its  derivative  u(t)  are  uniformly  bounded. 

In  the  subsequent  sections  It  will  be  shown  that  the  composite 
state/parameter  error  vector  Z(t)  of  several  adaptive  schemes  Is  governed  by 
the  same  general  equation 

Z(t)  -  A^(t)Z(t)  +  En(t)  +  F0(t)  +  Fe^^Q  (2.13 

where  0(t)  satisfies 

9(t)  -  A0(t)  +  Rn(t),  0(0)  -  0,  (2.14 

T  •  T 

E"  (-H  •  0]  and  F,  Q,  R,  and  A  are  constant  matrices  and  ReX(A)  <0.  We 
point  out  that, although  the  homogeneous  part  of  (2.13)  Is  linear,  this  system 
Is  not  Input  to  state  linear  because  A  (t)  depends  on  x  which  In  turn  depends 
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on  n.  For  this  reason  we  Introduce  the  notion  of  dominantly  rich  Inputs  in 
Section  2.3.  We  then  use  this  notation  to  establish  u.a.s.  of  the  homogeneous 
part  of  (2.13)  for  each  particular  adaptive  scheme.  Under  this  condition  the 
following  lemma  furnishes  a  bound  for  the  state/parameter  composite  error  Z(t). 
Lamma  2.1.1s  If  the  homogeneous  part  of  (2.13)  is  u.a.s.  then  Z(t)  is  bounded 


where 


lim  suplZ(t)ll  <  UY  —  ag(l+f) 
t^  ®2 


Y  ■  aup|u| 


(2.15) 


(2.16) 


.-l...-l_ .  IFJBRO  1 


(2.17) 


and  mj^,  m2,  a,  f^^,  and  f2  are  positive  constants  obtained  from  the  state 
transition  matrices  of  (2.13),  and  A  as  follows 


-m  (t-T)  Aj(t-T)  __g(t-T)  .  A(t-T).  .  . 


-U.-VU-1/  -flft-T)  A  (t-T) 

«♦  (t,T)l  <  m,e  ^  '  le  ^  I  <  ae  ®  <  f.e 

n  1  i 


Proof.  From  (2.13)  and  (2.18)  we  have 


-m.t  t  -m- (t-t) 

IIZ(t)l  <  m  e  ^  IZ(0)l  +  /  m,e  ^  |Elln(T)|dT 
^  0 


(2.18) 


t  -m,(t-T)  t  -m-(t-T)  -f-T 

+  /  m,e  ^  IFH9(T)|dT  +  /  m  e  iFlf  e  ^iQldx.  (2.19) 

0  ^  0 


,(t-T)  -f-T 

^  I  ^  ^  I 


In  view  of  (2.8),  (2.10),  (2.14),  and  (2.18)  we  can  write 


t  _ 

ln(t)J  <  ae'®^^''an(0)l  +  /  ae  “  lA’^jYdt 

0 


(2.20) 


|Rlln(T)ldT 


(2.21) 


lie(t)l  <  f,  /  e 
^  0 


where 


Y  ■  sup  I u 1 ,  6  ■  min | X (A-) |  i  I A-^l  ,  f ^  I ^ (A) | . 

t  ^ 


■  »  *2 

Substituting  for  nCt)  and  9(t)  In  (2.19)  we  obtain 


(2.22) 


lZ(t)l  <  p^e  “2*"  +  +  P3e  +  UY  ^  f 

(IIHO  +  IFIIRI)  (2.23) 

where  p^^,  °3  independent  of  t  and  bounded  for  all  uS  [0,1].  Thus 

as  t-^*,  (2.15)  follows. 

It  is  useful  to  emphasize  the  significance  of  the  terms  appearing 
on  the  bound  (2.15) 

g(l+f)  ~  characterizes  the  gain  in  the  disturbance  path 

a,m^  -  characterize  the  decaying  properties  of  the  parasitic  modes  and 
the  modes  of  the  homogeneous  part  of  (2.13),  respectively 
m^  ~  the  convergence  rate  of  the  homogeneous  part  of  (2.13) 

Y  -  a  measure  of  the  frequency  content  of  the  input  signal 
u  -  the  speed  ratio  of  slow  vs.  fast  phenomena. 

In  the  subsequent  sections  UY  will  emerge  as  the  most  Important 
factor.  It  is  already  clear  from  (2.15)  that  if  Lemma  2.1.1  holds  and  y<0(-^) 
then  the  scheme  is  robust  in  the  sense  that  the  composite  error  tends  to  zero 
as  the  frequency  of  parasltics  tends  to  infinity,  that  Is  as  y-^^O.  However, 
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if  Y>0(^)  this  robustness  property  cannot  be  ascertained.  Further  discxission 
of  the  bound  (2.15)  and  an  example  are  given  in  Section  2.5. 


2.  2.  Dominantly  Hich  Inputs 

A  well  known  requirement  for  asymptotic  identification  is  that  the 
input  u(t)  be  "rich"  enough  to  "persistently  excite"  all  the  modes  of  the 
plant  [27-29].  This  is  required  when  all  the  modes  are  to  be  identified. 
However,  our  goal  is  to  identify  only  the  dominant  and  to  disregard  the 
parasitic  modes.  In  this  sitiiation  inputs  exciting  only  the  dominant  modes 

seem  more  appropriate .  We  therefore  introduce  the  notion  of  "dominantly 
rich"  Inputs  and  in  the  subsequent  sections  we  show  that  they  guarantee 
the  robustness  of  several  adaptive  schemes.  Although  our  analysis  is 
applicable  to  a  wider  class  of  signals,  such  as  those  defined  in  [29],  for 
clarity  this  presentation  is  restricted  to  inputs  of  the  form 

k 

u(t)  ■  Z  a,  sinu). t  (2.24) 

i-1  ^  ^ 

where  are  all  distinct  and  positive  and  a^y^O. 

Definition  2.2.1;  Signal  (2.24)  is  a  dominantly  rich  input  for  the  plant  (2.7), 
(2.8)  if  there  exists  u*>0  such  that  for  all  ue  [0,u*)  the  vector  [x^,u^]  ^ 
is  "persistently  spanning"  (ps)  in  the  sense  of  [29]. 

We  point  out  that  this  definition  differs  from  the  usual  sufficient 
richness  condition  for  (2.7),  (2.8)  which  requires  the  ps-property  of  the 
vector  [x^,x^,u^]^.  An  important  practical  question  is  whether  every  signal 
which  is  sufficiently  rich  for  the  nth  order  plant  without  parasitics  (2.12) 
is  also  dominantly  rich  for  the  plant  with  parasitics  (2.7),  (2.8).  As  the 


example  at  the  end  of  this  section  shows,  the  dominant  richness  can  be  lost 
If  some  frequencies  In  u(t)  are  >  OC-^) .  To  prevent  this  we  restrict  the 
frequencies  In  (2.24)  to  be  lower  than 

Lemma  2.2.1:  If  an  Input  (2.24)  with  a)^<0(-^)  i=l,...,k  is  sufficiently  rich 

for  the  plant  without  parasltlcs  (2.12)  then  It  Is  also  dominantly  rich  for 
the  plant  with  parasltlcs  (2. 7), (2. 8). 

Proof  ;  The  steady  state  response  of  x  to  the  almost  periodic  Input  (2.24) 

Is  almost  periodic,  that  Is 

k 

X  ■  E  c.  sinu^t  +  d.  cosoi.t.  (2.25) 

1-1  i  i  i  i 

The  dominant  richness  is  achieved  if  [x^,u^]  is  pa.  For  this  it  is 
sufficient  [29]  to  show  that  for  all  and 


px  +  qu«0  p«0,  q-0. 
Substitution  of  (2.24)  and  (2.25)  into  (2.26) 


k 

E 

1-1 


((pCj^  + qa^)sinu^  t  +  pd^^  cos  u>^t)  -  0, 


implies  that 


(2.26) 


(2.27) 


pc^  +  qa^  -  0,  pd^  -  0,  1-  1,2 . k. 


(2.28) 


Following  the  same  procedure  as  in  [29]  we  use  (2.28)  and  the  relation 

Cl  +  jdi  -  (ju^I-A)'’^(B  +  uH(yju>^I-Aj)"^A’^B^ja)^)a^  (2.29) 

to  obtain 


p(a>^I  +  A^)"^(B-yAH(u^ojJl  +  Aj)"Hj+U^<i»jH(u^a)Jl  +  A^)"^A"^j)  -  0  (2.30) 


which,  after  grouping  the  0(y)  terms  becomes 
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p(«Jl+(A^+0(u))^)"^(B^  +  0(u)+0(w^u)J))  -0.  (2.30) 

By  Lemna  4  in  [29] ,  If  the  pair  (A  *l'0(u),  B  +  0(ii)  +  OCy^u^))  is  controllable 

o  o  1 

then  (2.31)  implies  p^O.  As  long  as  u .  <  0(— )  the  controllability  of  (A  ,B  ) 

X  U  o  o 

guarantees  that  there  exists  u*>0  such  that  for  all  ye  [0,y*)  the  pair 

2  2 

(A  +0(y),  B  +0(y)+0(y  u.))  is  controllable  [30^1].  Then  p«0  and  hence 
o  o  1 

q*0  in  view  of  (2.28). 

Lemma  2.2.1  excludes  input  frequencies  of  0(-^)  and  higher.  It  is  of 
interest  to  Illustrate  that  in  the  presence  of  parasltlcs  high  input  fre¬ 
quencies  can  Indeed  destroy  the  dominant  richness  of  the  Input.  Consider 


X  -  -a^x  +  bj^u  +  aj^2^  (2.32) 

yn  ■  -a^n  +  yb2U  (2.33) 


where  a^,a2»u)>0.  The  input  u^sinut  is  sufficiently  rich  for  the  dominant 
part  (2.32)  without  parasitics  (n*0).  For  this  input  to  be  dominantly  rich 


the  eqxiatlon 


P 


1 

2  2 
(o)  +a‘) 


^*12^2^2 

,2  2^  2 
y  u  +8^ 


2  . 

“  ^12°2 

^  2  2  .  2. 

(y  u 


(2.34) 


should  imply  p  *  0 .  However  this  is  not  so  since  a  frequency  u  zeroing  the 

expression  in  the  brackets  can  be  found.  For  example,  if  a^*l, 

a2“l,  b2"-l,  bj^“0.5  then  given  any  ye  [0,y*)  there  exists  an  u>0(~)  such 
2  2 

that  y  u)  *  l+2y  and  the  term  in  the  brackets  is  zero.  For  this  reason  the 
input  frequencies  are  restricted  to  u<0(~). 

Let  us  now  show  that  the  dominant  richness  of  u  assures  the 
robustness  of  the  simple  identification  schera  [2], [29].  The  problem  Is  to 
identify  the  dominant  part  of  the  plant,  that  is  the  pair  (A,B)  in  (2.7), 


when  the  dominant  state  x  Is  available  for  measurement.  We  disregard  the 
presence  of  the  parasitic  Input  Hn  In  the  design  of  the  Identification 
algorithm,  that  Is  we  assume  the  nth  order  model 

X  ■  K(x  -x)  +  A  (t)x  +  B  (t)u  (2 

in  m  m  n 

where  K  Is  a  stable  matrix.  As  In  [2]  the  adaptive  laws  for  adjusting 


A  (t)  and  B  (t)  are 
m  m 


^  =  -rexT 


'9  -  -Feu^  U 

where  ■  A  (t)-A,  'y -  B  (t)-B,  and  e*x  -x  are  the  parameter  and  state 
in  n  in 

T 

errors,  respectively,  and  r»r  >0.  When  the  adaptive  algorithm  (2.35)  to 
(2.37)  Is  applied  to  the  actual  plant  (2.7),  (2.8),  then  the  following 
set  of  equations  describe  the  behavior  of  the  Identification  scheme  In  the 


presence  of  parasltlcs 


e  ■  Ke  +  <^x  +  yu  -  Hn 


«  -Fex^ 


y  -  -Feu^. 


Defining  Z(t)  -  [e'^  ,x'^  ,t'^1  where  x*  . ^'*'l’'^2 . 

,  'fj  are  the  1th  rows  of  and  V,  respectively,  and  denoting 


[I 
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[="l‘ 

- 1 

M 

to 

•  F  - 

uY2 

• 

• 

u 

• 

• 

• 

• 

UY 

n 

L  J 

(2.42) 


where  Yj  is  the  jth  row  of  F,  equations  (2.38)  to  (2.40)  become  (2.13)  with  F=0. 
To  analyze  the  stability  of  the  Identification  scheme  in  the  presence  of 
parasltlcs  we  first  establish  the  u.a.s.  of 


Z(t)  -  A  (t)Z(t)  (2.43) 

n 

the  homogeneous  part  of  (2.13).  Although  (2.43)  is  formally  similar  to  the 
error  equation  for  the  nth  order  system  without  parasltlcs,  the  crucial 
difference  is  that  A^(t)  in  (2.43)  depends  on  n.  For  this  reason  the  u.a.s. 
of  (2.43)  does  not  immediately  follow  from  [2],  [29]  and  is  established  here. 
Theorem  2*2.1;  If  u(t)  is  dominantly  rich  for  the  system  (2.7),  (2.8)  then 
(2.43)  is  u.a.s.  and  by  Lemma  2.1.1  the  composite  state /parameter  error  Z(t)  is 

bounded  by  ^ 

lim  suplZ(t)l  <  PY  —  cig.  (2.44) 

®2 


Proof ;  Define 


- 1 

O 

’m  o' 

P  - 

>  0  ,  Q  - 

0  I 

0  0 

tm 

where  M ■ -(K^F  +  FK)  >  0.  Since  the  derivative  of 

V(Z)  -  Z^PZ 


(2.45) 


(2.46) 


Me 


for  (2.43)  is 


V(Z)  -  -Z^QZ 


(2.47) 


system  (2.43)  Is  uniformly  stable.  From  (2.47)  and  (2<43)  we  conclude  that 
e-+-0,  Z-^0  as  t-*-«  and  (2.43)  reduces  to 

+  y^^u  ■  0,  i»l,2,...,n  (2.48) 

where  i,  and  are  the  corresponding  values  of  e,  and  4'^ 
associated  with  the  homogeneous  part  (2.43). 

In  Lemma  2.2.1  we  have  established  that  when  u(t)  is  dominantly  rich, 

then  [x  ,u  ]  is  ps.  Thus  (2.48)  implies  (^^“0,  i*l,2 . n  which 

proves  that  (2.48)  is  u.a.s.  Once  u.a.s.  of  the  homogeneous  part  (2.43)  is 
established  Lemma  2.1.1  furnishes  the  bound  (2.44). 

In  conclusion,  the  identification  of  the  dominant  part,  which  is 
a  type  of  reduced  order  identification,  differs  from  the  full  order  identifica¬ 
tion  [2], [29]  in  the  dependence  of  x  and  hence  A  (t)  on  n  and  the  presence 

n 

of  the  forcing  term  Hn  in  (2.38). 

Furthermore,  instead  of  the  full  state,  only  the  dominant  state 
appears  in  the  steady  state  equation  (2.48).  '  The  dominant  richness  condition, 
which  avoids  excessive  excitation  of  the  neglected  fast  modes,  guarantees  the 
u.a.s.  However,  the  asymptotic  identification  is  not  achieved  since  the 
forcing  term  Hn  produces  a  steady-state  parameter  error.  The  identification 
is  robust  in  the  sense  that  the  parameter  error  is  0(p)  and  tends  to 
zero  as  u->-0,  that  is  as  parasitics  disappear. 
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2.3.  Reduced  Order  Adaptive  Observers 

If  instead  of  the  dominant  state  x  only  the  output  Is  available  for 
measurement,  the  problem  of  Identifying  the  parameters  and  estimating  the 
dominant  state  x  for  the  plant  (2.7),  (2.8)  simultaneously  Is  a  reduced 
order  adaptive  observer  problem.  Here  we  assume  that  the  output  Is  of  the 
form 

y  -  c'^^x  (2.49) 

that  is  the  parasltlcs  are  only  weakly  observable  from  (2.49)  since  the 
dependence  of  y  on  the  parasitic  modes  is  0(y).  The  case  of  strongly  observ¬ 
able  parasltlcs,  requires  a  different  type  of  analysis  and  is 

postponed  until  Section  4. 

We  assume  that  several  different  nth  order  adaptive  observers 
[32-34]  for  SISO  systems  have  been  designed  to  identify  the  triple  (A,B,C) 
and  estimate  the  dominant  state  x  of  (2.7),  (2.8)  by  disregarding  the 
parasltlcs.  Then  we  Investigate  the  behavior  of  these  schemes  when  applied  to 
the  actual  plants  with  parasltlcs. 

We  start  with  two  different  types  of  Minimal  Form  Adaptive 
Observers  [32], [33].  We  have  shown  in  Appendix  A  that  their  stability 
properties  in  the  presence  of  parasltlcs  are  determined  by 

e  -  Ke  +  d((jFv+ 4’^q)  -  Hn  (2.50) 


“  ®1 


(2.51) 


where  'P  are  the  parameter  Identification  errors,  e^  is  the  output  error, 
e  is  a  function  of  the  observation  error  and  v,q  are  the  filtered  output  and 
input  signals  respectively.  To  apply  Lemma  2.1.1  we  rewrite  (2.50)  to 

(2.53)  in  the  form  of  (2.13)  with  F-0  by  defining 


Z(t)  -  [e’^,4.'^,'?'*^]^ 


- 1 

H 

dqT 

i 

“ 

A^(t)  ■ 

-rv  1 

I  0 

-Mq  1 

(2.54) 


(2.55) 


Then  the  uniform  stability  of  the  homogeneous  part  of  (2.13)  can  be  shown  by 
using  the  same  Lyapunov  function  as  in  the  case  without  parasltlcs  [32], [33] 
and  it  is  not  Influenced  by  the  fact  that  A^(t)  depends  on  v,  which  in  turn 
depends  on  n  through  x^.  A  sufficient  condition  for  u.a.s.  of  the  homogeneous 
part  of  (2.13)  is  that  [v  ,q  ]  be  ps.  As  in  Lemma  2.2.1,  if  u(t)  is  dominantly 

rich  then  [v^,q^]  is  ps  and  therefore  the  homogeneous  part  of  (2.13)  is 

u.a.s.  This  result  is  stated  as  follows. 

Corollary  2.3.1;  If  u  is  dominantly  rich  then  the  homogeneous  part  of  (2.13) 
is  u.a.s.  Furthermore  the  composite  error  vector  Z(t)  of  (2.13)  is  bounded  by 

“l 

lim  suplZ(t)l  <  UY  “  og.  (2.56) 

t-»-  "  “2 

Next  we  consider  a  Non-minimal  Adaptive  Observer  [2], [33]  whose 
behavior  in  the  presence  of  parasltlcs  is  analyzed  In  Appendix  B  and  is 


described  by 
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“^1®1  d«'^v+4^q) -hF0-h^e^^x(O) 

(2.57) 

-re^^v 

(2.58) 

-Mej^q 

(2.59) 

where  e^  is  the  output  error,  (^,  are  the  parameter  identification  error, 

V  and  q  are  functions  of  the  states  and  the  input  of  the  non-mlnlmal  repre¬ 
sentation  of  the  plant  and  6  satisfies  (2.14) .  Defining 

Z(t)  -  (ej^,*T,«{rT]T  (2.60) 

A^(C)  - 


1  dv^  1  dq"^ 

_  —  —  ~  1—  —  _ 

•-h^- 

■x(0) 

-Tv 

1 

1 

1  0 

,  F  - 

.  Q  - 

-Mq 

1 

0 

0  , 

(2.61) 


we  rewrite  (2.57)  to  (2.59)  in  the  form  (2.13)  with  E“0  and,  similarly,  as 
before,  obtain  the  following  result. 

Corollary  2.3.2:  If  the  input  signal  is  dominantly  rich  the  homogeneous  part  of 
(2.13)  is  u.a.s.  and  the  composite  error  Z(t)  in  (2.13)  is  bounded  by 

“l 

lim  sup8Z(t)ll  <  yy  —  age  (2.62) 

t->»  ®2  * 


Finally  we  analyze  the  Parameterized  Observer  [34]  whose  state 
observation  error  e  and  parameter  error  4p  is  shown  in  Appendix  C  to  satisfy 


e  -  -M(t)GM^(t)CC^e  +  (AM(t)  +  (Iy,Iu))Ap  -  Ae^^e(O)  -  A9  -  Hn  (2.63) 

Ap  ■  -GM^(t)CC^e. 


(2.64) 


Combining  (2.63)  and  (2.64)  the  equation  for  the  composite  error 
Z(t)  ■  [e^.Ap^]^  Is  of  the  form  of  (2.13)  where  [-A^  0^],  [e^(0),0^], 

R  *  H  and 


*n<'> 


-M(t)GM'’^(t)Cc’^ 

T  T 
-GM  (t)CC 


AM(t) +  [Iy,Iu] 
0 


(2.65) 


In  this  case  the  dominant  richness  appears  Indirectly  as  a  sufficient 
condition  for 

t+T 

0<kj^I</  M’^(T)CC%(T)dT  <  k^I  forallt>0.  (2.66) 

Constants  k^,  k^,  and  T  rate  the  degree  of  linear  Independence  of 
the  components  of  the  vector  M^(t)C.  In  the  absence  of  parasitics  a 
sufficient  condition  for  existence  of  k^»  k2>  end  T  Is  that  u  Is  a  sum  of 
sinusoids  with  at  least  n  distinct  frequencies  [34].  Then  [x  ,u]  Is  ps  and 
therefore  the  components  of  M^(t)C  are  linearly  Independent.  In  the  presence 
of  parasitics  this  nth  order  richness  of  u  and  hence  the  linear  Independence 
of  M  (t)C  can  be  destroyed.  To  prevent  this  we  assume  that  the  Input  Is 
dominantly  rich. 

Theorem  2.3.1:  If  the  input  u  is  dominantly  rich  then  the  homogeneous  part  of 
(2.13)  Is  u.a.s.  and  Z  Is  bounded  by 


lim  suplZI  <  UY  — ^  ag(l+f)  (2.67) 

t-K.  ®2 


f 


kj^min  X[G] 

(1  +  nk,  max  X[G]) 
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Proof !  Using  the  same  proof  as  In  Theorem  1  of  [34]  we  can  show  that  if  there 
exists  constants  k2>  and  T  such  that  (2.66)  Is  satisfied  then  the  homo¬ 
geneous  part  of  (2.13)  Is  u.a.s.  with  a  rate  of  convergence  no  less  than  m^ 
where 

k.  min  X  [6] 

m-  - - - - b*  (2.68) 

(l  +  nk^maxXlG])'^ 

Once  the  u.a.s.  of  the  homogeneous  part  of  (2.13)  is  established 
Lemma  2.1.1  furnishes  the  bound  (2.67). 


2.4.  Discussion  and  Example 

Adaptive  schemes  designed  for  plants  without  parasltlcs,  but  applied 
to  plants  with  parasltlcs  are  considered  robust  If  the  error  due  to  parasltlcs 
Is  0(y).  This  means  that,  first,  the  error  Is  bounded,  and  second.  It  tends 
to  zero  as  parasltlcs  vanish,  l.e.  as  u-^0.  Among  the  factors  common  to 
all  the  error  bounds  derived  in  this  paper, 

®1 

UY-^og,  (2.69) 

®2 

factor  UY  determines  the  robustness  property.  If,  by  a  choice  of  high 
frequency  input,  Y*sup|u|  Is  made  0('^)  or  larger,  then  liY^O(l)  and  the 
error  bounds  do  not  tend  to  zero  as  u-^-O.  This  will  be  the  case  If  u(t) 
contains  terms  sin  -^t  where  a  Is  a  constant,  that  Is  If  Input  frequencies 
are  In  the  range  of  parasitic  frequencies.  For  such  Inputs  (2.69)  does  not 


exclude  the  possibility  that  the  state/paraineter  errors  may  be  significant. 

In  fact  If  u  contains  frequencies  of  0(~)  the  fast  state  n  will  be  of  0(1), 

Hr)  will  be  persistent  for  all  u  and  therefore  robustness  will  be  lost.  Our 
results  also  show  that  the  presence  of  parasltlcs  may  cause  loss  of  richness. 
If  the  richness  Is  achieved  with  frequencies  In  the  parasitic  range.  On  the 
other  hand  the  dominantly  rich  Inputs  guarantee  the  robustness  property. 

Let  us  Illustrate  this  and  other  aspects  of  the  derived  error  bounds  using  the 


plant 


r  -5  ll  [".91  ri.45‘ 

X  ■  X  +  “ 

.-10  Oj  L.sJ  12-25  _ 


WXg  ■  -4Xj  -  2u 
y  -  [l,0]x 

which,  upon  the  transformation  r)~x^  +  0.5u  becomes 


(2.70) 


(2.71) 


(2.72) 


-5  ll  fl]  [0.5 

X  ■  X  +  u  + 

.-10  oJ  [2] 


vif)  ■  -An  +  u0.5u 


y  -  (l,0]x. 


(2.73) 


(2.74) 


(2.75) 


Note  that  this  Is  an  example  of  the  plant  (2.1),  (2.2),  (2.49)  with  ^21*^ 
hence  L“0,  A22"A^,  Suppose  that  we  neglect  the  parasltlcs  and 

design  a  minimal  form  adaptive  observer  [32]  for  the  dominant  plant  (2.73), 


(2.75)  with  n“0,  that  is 


-6  1 


-8  0 


8  -  a.,(t) 


6-a  (t)]  b.(t) 

x+  y  +  u-e 


b,(t) 


^  -420v,v,+75v? 


r  ° 

:  1-A20v,v, 


,+75v 


I 


d 

■ 


(• 
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a2(t)  ■  75v2ej^,  bj^(t)  -  -Sqj^e^,  b2(t)  -  -7.8q2ej^ 

(2.77) 

V2  ■  -3v2  ^2  "  ”^*^2  *  ‘li  (2.78) 

A  A  ^ 

where  x  Is  an  estimate  of  x  and  a^(t)»  b^(t)  1*1,2  are  the  estimates  of  the 

unknown  parameters.  The  adaptive  observer  (2.76)  to  (2.78)  is  now  applied  to 
the  actual  system  (2.73)  to  (2.75)  with  parasitics  and  (2.73)  to  (2.78)  are 
simulated  on  a  digital  computer. 

The  dependence  of  the  error  bound  on  u  Is  Illustrated  in  Figs.  2.2 

and  2.3,  For  u  =  0.2  and  input  u  =  5  sin  2.5t  that  Is  y  *  18.2,  the  observation 
error  e2  is  relatively  small.  However,  the  parameter  errors  are  significant: 
10.4%  for  aj^(t)  and  12%  for  b^(t).  Reduction  of  u  by  a  factor  4,  that  is 
U*0.0S,  results  In  a  reduction  of  the  parameter  errors  by  approximately  the 
same  factor  as  shown  in  Fig.  2.3b,c.  The  observation  error  is  almost  zero 
in  this  case  (Fig.  2.3a).  Thus  for  a  fixed  frequency  the  parameter  error  is 
0(vi)  as  predicted  by  (2.69). 

To  examine  the  effect  of  Y*sup|u|  on  the  error  bound  the  value  of 
p  Is  kept  the  same  as  in  Fig.  2.3,  but  the  amplitude  of  the  higher  frequency 
Is  Increased,  u*  5  sin  t  + 15  sin  2.5t,  that  is  y  is  Increased  to  y*42.5.  As 
shown  in  Fig.  2.4  increasing  y  by  a  factor  of  2.3  results  in  an  increase  of  the 
parameter  error  by  a  factor  of  about  10.  Moreover,  the  observation  error, 
although  bounded,  is  oscillatory  and  not  close  to  zero. 

An  even  more  critical  way  to  change  the  frequency  content  of  the 
input  signal  is  to  increase  the  higher  frequency  ou*2.5  tenfold  to  u*25, 

u*  5  sin  t -f  S  sin  2St.  Decreasing  the  value  of  u  by  the  same  factor,  from 
u*0.2  to  u*0.02,  does  not  reduce  the  identification  error  which  in  Fig.  2.5 


is  al^oue  3  times  higher  chan  in  Fig.  2.2.  This  illustrates  that  input  fre¬ 


quencies  of  0(-^)  may  cause  state/parameter  errors  which  persist  or  even 


increase  as  u-»-0.  Hence  loss  of  robustness  and  a  possibility  of 
instability. 

A  trade  off  in  selecting  the  input  is  apparent  from  this  discussion. 
High  frequencies  are  to  be  avoided  since  they  excite  the  parasitic  modes .  On 
the  other  hand  low  frequencies  affect  the  convergence  rate  adversely.  It  is 
appropriate  to  use  a  dominantly  rich  input  with  a  low  value  of  y  still  giving 
an  acceptable  race  of  convergence. 


The  improvement  of  the  bound  by  making  the  convergence  rate  m2  larger 


may  require  trial  and  error  selection  of  the  adaptive  gains .  An  exception  is 
the  parameterized  adaptive  observer  where  the  expression 


m 


kj^min  X[G] 


^  (1+ nk2  max  X[G])^ 


(2.79) 


gives  more  information  about  the  dependence  of  the  error  bound  on  other 
quantities.  In  this  case  m2  is  the  largest  when  maxX[G]“  minX[G].  Note 
from  (2.79)  that  the  error  bound  will  be  higher  when  the  order- n  of  the 
dominant  part  of  Che  plant  is  higher. 
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3.  DISCRETE-TIME  PARALLEL  AND  SERIES  PARALLEL  IDENTIFIERS 

3.1.  Characterization  of  the  Mismatch 

A  class  of  discrete  systems  possessing  a  two-time  scale  property  has 
the  form  [35] 


x(k+l) 

-  Aj^j^x(k)  +  pAj^2*W  +  Bj^u(k) 

(3.1) 

z(k+p 

-  A2^x(k)  +  MA222(k)  +  B2U(k) 

(3-2) 

y(k) 

-  Cj^x(k)  +  y.C2Z(k) 

(3.3) 

where  exists  and  p,  Is  a  small  positive  parameter.  The  states  x(k)  and  z(k) 
are  n  and  m  vectors  respectively  and  u(k)  and  y(k)  are  the  scalar  Input  and 
output  respectively.  The  restriction  of  the  output  of  the  model  to  be  of  the 
form  (3.3)  makes  the  parasltics  weakly  observable  in  the  plant  output.  This 
weak  observability  assumption  is  found  to  be  crucial  in  estabilishing  the 
robustness  of  the  adaptive  schemes,  i.e.  as  (modeling  error  -+•0)  the 

composite  identification  error  ^0.  The  case  of  strongly  observable  parasitics 
requires  a  different  kind  of  treatment  and  it  will  be  discussed  in  Section  4. 


We  use  the  tranformation 

n(k)  ■  z(k)  +  Px(k)  -  B£u(k)  (3.4) 

to  obtain  the  representation  (3.5)  tr  (3.7) 

x(k+l)  -  A  x(k)  +  B  u(k)  +  pHn(k)  (3.5) 

S  8 

ri(k+l)  ■  MA^n(k)  +  Bj(u(k)-u(k+l)  ]  (3.6) 

y(k)  ■  CgX(k)  +  uC2B£u(k)  +  yC2n(k)  (3.7) 


where  P  satisfies 
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^21  ■  ^^^22^  “  ^***^12^  “  ° 

«ad  H  -  A^2.  A3  -  B3  -  +  p.A^2®f*  '  h  '  ^  *  ^22 

PAj^2»  Bf  “  (I  “  y^Aj)  ^(B2  +  PBj^).  Representation  (3.5)  to  (3.7)  is  found,  to 
be  convenient  for  characterizing  the  error  bounds  of  the  adaptive  schemes 
considered.  Without  loss  of  generality  let  us  assume  that  the  model  of  the 
dominant  part  of  the  plant  is  in  the  observable  canonical  form,  i.e. 


.  -  [1  0  ....  0] 


and  obtain  the  following  representation  for  the  plant 


X  (k+1)  ■  A  X  (k)  +  b  u(k)  +  wH  n(k) 
s  s  s  s  s 

(3.9) 

n(k+l)  -  yA^n(k)  +  B^[u(k)-u(k+l) ] 

(3.10) 

y(k)  -  (1  0  ...  0]Xg(k).-  C3X^(k) 

(3.11) 

also  shown  in  Fig.  3.1 

where  b^  -  Bg-uagC2Bj  +  u[C2Bj,0. .  .of  ,  -  H-agC2+y  [C2A^,0.  .  .0] 

Due  to  this  observable  canonical  form  we  can  write  (3.9)  to  (3.11)  as 


y(k)  -  Z  a.y(k-i)  +  Z  b.u(k-i)  +  y  (k) 
i-1  ^  i-1  i  « 

where  y^^(k)  is  the  unmodeled  part  of  the  plant  and  is  given  by 


(3.12) 


y^(k)  -  y  Z  h^n(k-i). 

X  X 


(3.13) 


and  h.  is  the  ith  row  of  H  . 
1  s 


Assume  now  that  an  Identifier  is  designed  for  the  nth-order  dominant  part 

T  T 

to  identify  the  vectors  a  ■  [a,  ...  a  ]  ,  b  ■  [b.  ...  b  ]  under  the 
assumption  that  there  are  no  parasltlcs,  l.e.  n(k)>0.  This  Identifier  is 
then  applied  to  the  actual  system  with  parasltlcs.  Our  purpose  Is  to  examine 
the  robustness  of  "parallel”  and  "series  parallel"  Identifiers  In  the  presence 

of  the  parasitic  Input,  obtain  bounds  for  their  parameter  and  output  errors, 

and  ther^y  examine  their  performances.  Since  the  unknown  parameters  s  ,b 

s  s 

are  functions  of  |J.  to  each  value  of  p.  there  corresponds  a  different  parame¬ 
terization. 

The  following  assumptions  are  made  throughout  this  section: 

(Al)  A  and  A,  are  stable, 
s  f 

(A2)  The  order  n  of  the  dominant  part  of  the  plant  is  known. 

(A3)  The  triple  (A  ,b  ,C  )  Is  both  completely  controllable  and  completely 

s  s  s 

observable. 

(A4)  The  only  available  signals  are  the  Input  uCO  and  output  y(>). 

(A5)  u(k)  and  u(k)-u(k+l)  (and  therefore  y(k)  due  to  (Al))  are  bounded  for 
all  k. 


3.2.  Reduced  Order  Parallel  Adaptive  Identifier 

Equations  for  an  nth-order  parallel  identifier  [6]  are 
summarized  below.  The  estimation  model  (or  parallel  system)  is 
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y®(k)  »  Pp(k-l)0p(k-l)  (3.15) 

where 

p]J(k)  -  (a  (k)  ...  a  (k)  b,(k)  ...  b  (k)],  (3.16) 

pi  n.  1  n 

ep(k-l)  -  [yp(k-l)  ...  yp(k-n)  u(k-l)  ...  u(k-n)]^,  (3.17) 

and  y^(k),  yp(k)  are  the  a  posteriori  and  a  priori  output  of  the  estimation 
model,  respectively,  at  the  Instant  k.  The  Integral  adaptation  algorithm  is 
[61 

F^(k-I)e^(k-I)v°(k) 

"  P-(k-l)  + - ^ 2 2 -  (3.18) 

^  P  l+6i(k-l)F^(k-l)e  (k-1) 

P  P  P 

Vp(k)  -  y(k) -Pp(k-l)ap(k-l)  +  ^?^c^ty(k-i)-yp(k-i) ] ,  (3.19) 

where  Fp(k-l)  can  be  a  decreasing  positive  definite  matrix  given  by 


F^  (k-1)  9^  (k-1)  e;:  (k-1)  F  (k-1) 
F  (k)  -  F^(k-1)--E - — E - 2 - 2 - 

^  P  X+9^(k-l)F  (k-1) 9  (k-1) 

P  P  P 


X  >  0.5 


(3.20) 


where  F  (k  )  Is  sn  arbitrary  positive  definite  matrix,  r  The  constar.ts  c.  ,..,c 
p  o  ■  X  •  n 

are  preselected  such  that 


H(2)  - 


^  -i 


1-  Z  a. z 
i«l  i 


-i 


_1_ 

2X 


(3.21) 


is  strictly  positive  real,  i.e.  Re{H(z)}>0,  vizI  ■  1. 

Theorem  3.2.1:  With  u(k)  satisfying  a  persistently  exciting  condition  [36], 
the  adaptive  identifier  (3.14)  to  (3.21)  is  stable  in  the  presence  of  a 
parasitic  input  (3.13)  in  the  sense  that  the  composite  error 


35 


a-  [a.  a,  ...  a  1^,  b-[l  0  ...  0]^,  c-tc,  c.  ...  c  ]^,  (3.29) 

j.  z  n  1  z  n 


1  0  .  0 

•  • 

0  1  *•*.  : 

•  ••  *•.  : 

•  •  •  * . 

•  *.  •  •  * 

.  •  ♦  . 

0  ...  0  1  0 

Not.  th.t  (3.26)  1.  not  llo«r  ..  It  .pp«r.  .i„c.  A^(k)  depends 

on  e„(k)  which  Includes  y  (k)  which  is  in  Z  (k) .  As  in  [21]  define  A  (k)  to 
P  P  P  P 


be  Ap(k)  with  9p(k)  replaced  by 

0(k)  -  [y(k),y(k-l) . y(k-n+-l)  ,u(k-l)  , . . .  ,u(k-tt+l)  ]  (3.31) 

and  F  replaced  by  f  (k)  with  0  (k)  replaced  by  6(k)  and  set 
P  P  P 

(Z  (k)  ,k)  -  A  (k)  -  A  (k)  .  (3.32) 

r  H  It 

Note  that  A  (k)  is  independent  of  Z  (k)  and  that 
P  P 

9(k)  -  9p(k)  +  [ep(k)  0]'^.  (3.33) 

Then  (3.26)  can  be  rewritten  as 

Zp(k+1)  -  Ap(k)Zp(k)  +  iKZp(k),k)Zp(k)  +  Bp(k)y^(k+1) .  (3.34) 


Since  A  is  identical  to  F  in  [36]  with  A  -  «*  and  B  is  bounded,  then  if  u 
P  P 

is  sufficiently  rich  (in  order  to  insure  the  exponential  stability 
Zp(k+1)  -  Ap(k)Zp(k)  [36])  then 


* 

With  the  retention  of  positive  definiteness  and  invertlblllty  by 
Fp(k)  in  (3.20),  the  time-variations  of  Fp(k)  do  not  violate  the  uniform 
observability  requirements  of  [36]  nor,  given  the  global  asymptotic  stability 
proven  in  [6]  for  (3.20),  the  conclusion  of  exponential  convergence.  Note 
Fp(k)  does  not  retain  positive  definiteness  for  all  A,  e.g.  A>1  even  given 
sufficient  excitation  period. 


(3. 


lim  suplZ  (k)l  <  S(Z  (k  ),p)lim  sup  y  (k+1)  , 
P  P  o  u 


where  S(Z  (k  ),p)  is  a  boxinded  function  of  the  Initial  parameter  and  output 
P  o 


estimate  error  Z  (k  }  and  the  exponential  convergence  rate  p  In  the  absence 
P  o 


of  parasitic  Inputs. 

From  (3 . 10) 


k-1 


n(k)  ■  «>(k,k  )n(k  )  +  I  (^(k,l+l)B-[u(l)-u(l+l)  ] , 
o  o  f 


(3. 


where  ^(k^j)  Is  the  transition  matrix  of  (3. 10),  l.e.  (yA^)  .  Since  Is 

stable,  [37] 


l<fr(k,j)l  <  aj^(ya2)  for  •  >  >  0  and  1  >  02  >  0. 


Therefore,  In  (3.36), 


k-k 


k-1 


In(k)|  i  a.  (ya.)  ®lln(k  )|+  Z' a,  Cyo,)^”^"SB-|  lu(l)-u(l+l)  | .  (3. 


o  1-kj,  1  2 


Defining 


Y  ■  suplu(l)-u(l+l) 


converts  (3.38)  to 


k-k  o. 

In(k)l  <  a^(ya2)  °ln(k^)l  +iI^»BjIy. 


From  (3.13) 


V  (k+l)|  <  yn(  sup  lh.l)(  sup  ln(k-l+l)|). 


n>l>l  ^  n>l>l 


Defining 


h  -  sup  I h . 1 , 
n>l>l  ^ 


then,  since  0  <  <  1> 


11m  sup  I y  (k+1) [  <  ynyhl B-1 
k-*«»  “  ‘ 


1-ya, 


(3. 


(3. 


(3. 


(3. 


(3. 


(3. 


Substituting  (3.43)  Into  (3.35)  establishes  (3.23). 
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3.3.  Reduced-Order  Series-Parallel  Adaptive  Identifier 

Consider  now  the  series-parallel  (or  equation  error)  Identifier 
configuration  [6], [38].  The  parallel  adjustable  model  given  before  is 
replaced  by  a  series-parallel  adjustable  discrete  system 

y.(lc)  -  ?  i.(k)y(k-i)  +  Z  b,(k)u(k-l)  -  p^(k)0(k-l) 

y°(k)  -  p3(k-l)9(k-l) 
s  s 

where 

0(k-l)  «  [y(k-l),y(k-2),...,y(k-n),u(k-l),...,u(k-n)l^ 

and 

p3(k)  -  [a-(k),a,(k),...,£  (k),b,(k) . b^ (k)]. 

s  X  z  n  1  n 

The  adaption  algorithm  to  be  Investigated  Is  [6]* 

F^(k-l)0(k-l)v®(k) 

P  (k)  -  5  (k-1)  + - S- - S -  (3^48) 

®  ®  l+0^(k-l)F  (k-l)8(k-l) 

s 

v°(k)  «y(k)-y°(k)  (3.49) 

s  s 

where 

F  (k-l)9(k-l)9^(k-l)F  (k-1) 

F  (k)  -  F  (k-1)  - - = - § -  ,  X>0.5.  (3.50) 

X  +  0  (k-l)F  (k-1)  0  (k-1) 
s 

Theorem  3 .3.1;  For  the  series  parallel  identifier  of  (3.44)  to  (3.49)  the 
composite  error 

Z.(k)  -  [y(k)-y  (k),a.-a. (k),...,a  -a  (k),b,-b,(k) . b  -b  (k)]^  (3.51) 

s  sxi  nnii  nn 

is  hounded  for  all  k  and 

A 

For  X»1  the  algorithm  (3.48)  to  (3.50)  is  the  standard  least 
squares  algorithm  which  is  not  exponentially  stable  and  therefore  lags  the 
properties  proven  in  [36]. 


(3.44) 

(3.45) 

(3.46) 

(3.47) 


(3.52) 


11m 

It** 


suplZ  (k)l  <  pnyhlB 

8 


1  _ ± _  -2_t 

:  I  l-ya2j  [1-C2J’ 


where  &  Is  related  to  the  bounds  on  u,  y,  and  F  and  c.  and  c.  In  c. (c.) 
s  s  1  A  1.  z 

bound  the  k  to  k  transition  matrix  of  the  linear  homogeneous  Identifier  error 
o 

system  (l.e.  below). 

Proof:  The  Identification  error  system  obtained  In  Appendix  E  is 


where 


Z^(k+1)  -  A^(k)Z^(k)  +  B_(k)y  (k+1) 
s  s  s  s  u 


A,(k) 

s 


e'^(k)  J 


f  Fg(k)0(k)0'^(k)') 


l+0"-(k)F  0(k) 
s 


T 


I  - 


P.(k)0(k)0^(k) 

S 

l+0’^(k)F  (k)0(k) 
s 


B^Ck) 


_ 1 _ 

l+0^(k)F  (k)0(k) 
s 

-Fg(k)0(k) 

l+0'^(k)F  (k)0(k) 
s 


(3.53) 


(3.54) 


(3.55) 


The  exponential  stability  of  the  homogeneous  portion  of  (3.53)  is  proven  in 
[36]  (with  X =  “) .  Therefore,  the  bound  from  [21]  used  in  (3.35)  is  applicable. 
However,  since  in  this  case  the  y^^O  convergence  rate  of  IIZg(k)II  is 

determinable  directly  from  the  transition  matrix  of  (3.53)  and  is  not  a 
function  of  Zg(k^)  as  in  the  parallel  case.  Due  to  the  exponential  stability 
of  the  homogeneous  portion  of  (3.53) 


k-k 


IZ^(k)l  <  Cj^C2 


+^Z^c,c^"^"Sb  (DJlv  (i+l)«  . 
i-k-  12  s  u 


(3.56) 


From  (3.55),  due  to  the  boundedness  of  u,  the  stability  of  Ag>and  the 

boundedness  of  F  (k) ,  B  (k)  Is  bounded.  Define 
8  8 

6  ■  suplB  (1)1  <  •.  (3.57) 

®  1  * 

Using  (3,43)  and  (3.57)  In  (3.56)  results  In  (3.52)  as  k->-<». 


3.4.  Discussion  and  Example 

In  full-order  application  parallel  and  series-parallel  Identifiers 
differ  In  their  possible  unbiasedness  or  biasedness,  respectively.  In  the 
presence  of  zero-mean  output  measurement  noise  and  In  the  presence  or 
absence  of  a  priori  plant  Information  requirements,  respectively,  as  Indicated 
by  the  strict  positive  reality  requirement  In  (3.21)  for  parallel  Identifiers 
and  the  lack  of  any  such  requirement  for  series-parallel  Identifiers.  As 
shown  In  this  study  these  two  classes  also  demonstrate  differences  in 
reduced-order  application  since  the  bound  (3.23)  for  parallel  Identifiers  Is 
dependent  on  the  Initial  composite  error  and  the  "homogeneous"  algorithm 
convergence  rate  but  for  series-parallel  Identifiers  on  the  Identifier  error 
system  Input  matrix  norm  maximum  and  the  "homogeneous"  algorithm  convergence 
rate  but  not  the  Initial  composite  error.  However  both  bounds  share  the 
common  factor 

M>Y»8  (3.58) 

where  g  ■  h  I!  B-l|  _ ^ _  showing  that  both  identifiers  are  robust  with 

1  "  |J^2 

respect  to  unmodeled  parasltlcs  In  the  sense  that  as  |j,->0  the  error  bounds  go 
to  zero.  The  dependence  of  (3.58)  on  y  •  sup  |u(l-t>l)  -  u(l)  |  is  an  indicator 
of  the  Influence  of  the  input  characteristics  on  the  identification  results. 
The  bounds  Indicate  possibilities  for  reducing  the  error  by  a  proper  choice 
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of  the  input  signal.  That  is  the  identification  results  in  the  presence  pf 
parasitics  can  be  improved  by  choosing  an  input  signal  which  is  sufficiently 
rich  but  has  a  low  y  >  sup  |u(i  +  l)  -  u(i)|.  Some  of  the  terms  in  the  bounds 
obtained  are  difficult  or  even  impossible  to  calculate  since  they  depend  on 
the  characteristics  of  the  unknown  parasitics.  Thus  the  bound  can  only  be 
used  for  qualitative  analysis  and  not  as  a  practical  measure  for  evaluating 
identification  results  in  any  given  practical  situation.  Such  a  qualitative 
analysis  indicates  possible  Improvement  of  the  identification  results  by  a 
proper  choice  of  the  input  signal  and  an  increasing  rate  of  convergence  of  the 
algorithm. 

The  effect  of  p.,  the  input  characteristics  Y  and  the  initial  com¬ 
posite  error  on  the  identification  results  has  been  demonstrated  by  simulating 
the  following  example: 


x(k+  1) 

1 

2(k+  1) 

-0.2 

1 

B 

D 

jH 

B 

B 

0 

0 

0 

O.TjJ 

x(k) 


z(k) 


1 

0 

0.8 

0.1 


(3.59) 


y(k)  -  [1  0]  x(k) 


(3.60) 


which  can  be  expressed  in  the  form  of  (3.12)>  (3.13)  by  using  the  transforma¬ 
tion  (3.4)  i.e. 

y(k)  -  -0.2y(k  -  1)  +  0.48y(k  -  2)  +  (1  +  ^°'-Q’^3^)u(k  -  1) 


+  ( 


0.4u 

1-0.5^ 


+ 


0.  lu 

1  -  0.7p 


)u(k  -  2)  +  y^(k)  , 


(3.61) 


y,  (k)  -  p(n,  (k-  1)  +  O.Sn,  (k-  2)  +  no(k-  1)) 


(3.62 
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The  criteria  used  for  evaluation  of  the  Identification  results  are  the  para- 

“  ^2  '‘2  1/2 

metric  distance  defined  by  D  «  ^£^1^*1  ”  "*■  the 

output  error  e  ■  y(k)  -  y(k) .  For  the  Input  excitation  signal  we  have  chosen 
a  sequence  of  pulses  whose  frequency  and  amplitude  have  been  adjusted  to  give 
a  particular  value  of  Y< 

The  effect  of  the  Initial  error  on  the  performance  of  the  parallel 
Identifier  Is  shown  in  Fig.  3.2  and  Fig.  3.3.  That  is  by  increasing  the  norm  of 
the  Initial  error  from  1  to  IS  the  "convergent"  parametric  distance  Increases 
by  55.27.  (l.e.  from  0.087  to  0.135).  The  same  change  In  the  Initial  error 
has  no  effect  on  the  performance  of  the  series-parallel  Identifier  as  indicated 
by  Fig.  3.4  and  Fig.  3.5.  The  effect  of  |x  on  the  Identification  results  Is 
shown  In  Figs.  3.4  and  3.6.  By  increasing  v  by  5  times,  the  "convergent"  para¬ 
metric  distance  Increases  by  5.44  times  and  the  oscillations  of  the  output 
error  become  higher.  This  shows  that  when  the  parasltlcs  are  sufficiently 
fast  the  Identification  results  can  be  satisfactory.  The  effect  of  the  char¬ 
acteristics  of  the  Input  on  the  Identification  results  Is  deoonstrated  by  simu¬ 
lations  and  is.  shown  in  Figs.  3.4  and  3.7.  By  Increasing  y  *  sup|u(l+l)-u(l) | 
from  13  to  100  the  "convergent"  parametric  distance  Increases  by  5.47.  and  the 
output  error  deviations  from  zero  Increase  considerably. 

For  both  classes  of  algorithms  the  "homogeneous"  l.e..  In  the  absence 
of  parasltlcs,  convergence  rate  Is  a  complicated  function  of  Input  {u(k)}  com¬ 
position  and  the  step-size  matrix  F(k).  For  parallel  Identifiers  the  error 
smoothing  coefficients  c^  In  (3.21)  also  Impact  the  convergence  rate.  Clearly, 
further  refinement  of  reduced-order  usage  composite  error  bounds  will  require 
more  detailed  descriptions  of  the  convergence  rate  mechanics ,  which  is  an 
admittedly  [21]  difficult  task. 


the  Parallel  Identifie 
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Fig.  3.7  Identification  results  using  the  Series-parallel  identifier  for  a 
V  -  100  and  ||Z(Ko)||  -  1. 


0.1, 


4.  ROBUSTNESS  IMPROVEMENT  OF  IDENTIFIERS  AND  ADAPTIVE  OBSERVERS 


4.1.  Introduction 

In  Sections  2  and  3  It  was  shown  that  several  adaptive  schemes  are 
robust  with  respect  to  modeling  errors  consisting  of  fast  parasltlcs  which 
are  weakly  observable  In  the  plant  output.  In  this  section  we  show  that  the 
assumption  of  weak  observability  Is  crucial  for  robustness.  If  the  parasitic 
modes  are  strongly  observable  the  schemes  are  no  longer  robust.  That  Is  the 
error  due  to  model-plant  mismatch  can  be  large  even  when  the  fast  modes  are 
Infinitely  fast.  In  the  case  of  continuous-time  adaptive  schemes  we  show 
that  the  addition  of  a  low  pass  filter  at  the  output  makes  the  parasltlcs 
weakly  observable  and  hence  guarantees  the  robustness  of  the  enlarged 
scheme.  The  robustness  of  discrete-time  adaptive  schemes  with  respect  to 
strongly  observable  parasltlcs  can  be  guaranteed  by  modifying  the  adjustable 
model  and  the  adaptive  laws. 

4.2.  A  Model  with  Fast  Parasltlcs 

A  SISO  system  possessing  slow  and  strongly  observable  fast  parts  can  be 
represented  In  the  explicit  singular  perturbation  form  [25] 


^1*  ^2*f  ^l“ 

(4.1) 

A21X  +  A^jXf  +  b^u 

(4.2) 

(4.3) 

where  x  e  6  R®  and  p  Is  a  small  positive  parameter  associated  with 

the  presence  of  "parasitic"  elements,  such  as  the  time  constants,  masses,  etc. 


Without  altering  the  input-output  characteristics  of  the  system  we 


use  the  transformation 


n  ■  +  Lx  +  A^H^u 


to  obtain 


>>o"  *  ^2" 


pn  “  A-n  +  wb.u 


y  -  C^x  +  C2n-C2b^u 


(4.^ 


(4.5 

(4.6 

(4.7 


where  ,  A^-A22+  ^o  “ '*l“^l2^f^f  ’  *"o 


C1-C2I 


and  L  satisfies  the  algebraic  equation  [26] 


A22L  -  A2j^  +  "  '^^ll  “  °  (^  •  ‘ 

Assuming  that  the  pair  (C^,A^)  is  completely  observable  (4.5)  to  (4.7)  can  be 
written  in  the  following  "modal"  canonical  form 


(4.9 


no  ■  AjO  +  b^u 


(4.10 


Y  ■  CX  -  y,  C  »  [1,0,...,0]  (4.11 

where  A^*A^2“®gC2t  b*b^  +  a^C2bj,  h  “[11  ...  1],  0*^02 »  A  is  an 

(n-l)x(n-l)  diagonal  matrix  with  arbitrary  but  known  constant  and  negative 
diagonal  elements  — (1*2,..., n).  Here  the  fast  state  o  Is  defined  by* 


Note  that  for  the  initial  conditions  of  a  to  be  finite  as  iJ-+0  the 
initial  conditions  of  q  in  (4.12)  have  to  be  at  most  of  0(y). 


The  reduced  order  system  obtained  by  setting  ii"0  In  (4.9)  to  (4.11)  Is 


(4.13) 

(4.14) 


We  point  out  that  in  Section  2  we  asstimed  C2^0,  therefore  the 
parasltlcs  were  weakly  observable.  Here  we  allow  C^^Q.  This  means  that  as 
U-K)  the  transfer  function  of  the  system  becomes  proper.  Thus  the  input  depen¬ 
dent  signal  -C2b^u  is  present  in  the  dominant  part  of  the  plant  and  remains  in 
the  reduced  order  system.  Since  the  continue us -time  adaptive  schemes  are 
designed  for  strictly  proper  systems  the  appearance  of  -C2b^u  cannot  be  seen 
by  the  adaptive  algorithm.  The  purpose  of  this  section  is  to  investigate  the 
behavior  of  an  n-th  order  observer  with  respect  to  parasltlcs  when  €2)^0  and 
compare  with  that  in  the  case  of  €2*0.  We  will  first  show  that  the  appearance 
of  the  input  -C2b^u  in  the  reduced  order  system  is  equivalent  to  the  nonrobust¬ 
ness  of  the  adaptive  observer  with  respect  to  parasltlcs.  Then  we  will  show 
how  the  robustness  of  the  adaptive  scheme  can  be  reestablished  using  a  low-pass 
filter  at  the  output  of  the  plant.  The  following  assumptions  are  made: 


(Al)  A  and  A.  are  stable, 
o  t 

(A2)  The  order  n  of  the  dominant  part  of  the  plant  is  known. 

(A3)  The  triple  (A  ,b,C)  is  completely  controllable  and  observable. 

o 

(A4)  The  only  available  signals  are  u(t)  and  y(t). 

(AS)  u(t)  and  u(t)  are  piecewise  continuous  bounded  functions  of  time. 
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4.3.  Minimal  Form  Adaptive  Observer 

Suppose  that  an  n-th  order  adaptive  observer  Is  designed  to  estimate 
the  state  vector  X  of  the  dominant  part  of  the  plant  and  Identify  the  unknown 
vectors  a^  and  b.  It  Is  assumed  that  the  dominant  part  of  the  system  Is 
strictly  proper  and  the  presence  of  parasltlcs  Is  disregarded  during  the 
design.  The  stability  of  the  adaptive  observer  when  applied  to  the  real  plant 
with  a  proper  dominant  part  and  parasltlcs  Is  analyzed. 

The  algorithm  [32], [33]  for  the  n-th  order  adaptive  observer  based  on 
the  system  (4 .13) , (4.14)  without  the  signal  -C^b^u  Is  given  by  the  equations 
(4.15)  through  (4.16)  below.  The  observer  equation  is 

X  -  KX  +  [k-a(t)]y  +  fi(t)u  +  w  +  r  (4.15) 

y  -  (4.16) 

where  w  and  r  are  auxiliary  signals  formed  by  the  output  error  e^^  •  y-y, 
the  derivatives  of  the  parameter  error  components  and  the  components 


v^.y.  q^-u.  v^-  — y.  q^ 


of  the  vectors  v  and  q  as  follows 


(i-2 . u)  (4.17) 


*0  * 

fM 

w  ■ 

^2^ 

,  r  - 

cr 

CM 

< 

L  n  j 

“  j 

(4.18) 


Moreover 

strictly 


-X  I  h 

0  1--  • 

:  '  A 


Ti 


Is  stable. 


positive  real,  and  d*  (1  0 


T  -1  1 

the  transfer  function  C  (sI-K)  d*— Is 

s+X^ 

T 

...  0)  .  The  adaptive  laws  for  adjusting 


the  parameters  are 


where  it*a(t)-a  and  4'»b(t)-b. 
o 

It  can  be  shown  as  In  Section  2  that  the  behavior  of  the  adaptive 
observer  In  the  presence  of  parasltlcs  Is  equivalent  to  the  stability  of  the 
following  set  of  differential  equations 


£  ■  Ke  +  d[(J^v+f^q]-  (C2A^  + 

"l  -  ®1 
i  -  -FEj^v 

¥  -  -MEj^q 

which  using  the  composite  error  Z(t)  *  [e^ become 

Z(t)  ■  A  (t)Z(t)  +  Ha(t) 

n 


(4.21) 

(4.22) 

(4.23) 

(4.24) 


(4.25) 


where 


A^(t) 


I  dv' 


L-Mq  1 


0  I  0 


H  »  - 

0 

(4.26) 


The  following  theorem  establishes  the  boundedness  of  the  composite  error  Z(t) 
in  the  presence  of  parasltlcs. 

Theorem  4.  3.1:  If  u(t)  is  dominantly  rich  for  the  plant  (4.9)  to  (4.11)  then  the 

homogeneous  part  of  (4.25)  Is  u.a.s.  and  the  error  vector  Z(t)  Is  bounded  by 

m  a 

Urn  sup|Z(t)l  <  (IC,AJ  +wlA,,-a  C,l)  lb,l 

t-»-»  3  t  12  o  2  ®2  *2  ^ 

where  m.,  m.,  a.,  and  a.  are  positive  constants  in  the  bounds  |(^  (t,T)|  ^ 


(4.27) 
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-m-(t-T) 


^t-Tx 

-ao(T“) 


<  m^e  -  ,  l4)^(t,T)l  <  ^  ^  are  the  transition 

matrices  of  (4.25)  and  (4.10)  respectively,  and  Y“sup|u(t)|. 

Prop f ;  The  proof  of  Theorem  4.1.1  follows  directly  from  Lemma  2.1.1  and 
Theorem  2.3.1  of  Section  2  and  is  not  repeated  here. 

When  C2"0  the  algorithm  is  robust  with  respect  to  parasltics  because 
from  bound  (4.27)  the  identification  and  observation  error  goes  to  zero  as  u-^. 
However,  when  C2)^0  the  algorithm  is  no  longer  robust  with  respect  to  parasltics. 
That  is,  as  (l.e.  modeling  error  0)  the  composite  error, although  bounded, 
does  Qot  go  to  zero.  We  show  this  by  applying  the  adaptive  observer  to 
the  reduced  order  system  (4 .13) , (4 .14) .  In  this  case  the  stability  properties 
of  the  adaptive  observer  are  equivalent  to  the  stability  of  the  following  set 


of  error  equations 


E  -  Ke  +  dE^^v+i’iil  +C^b^u 

®1  “  ®1  " 
i  -  -FEj^v 
¥  -  -M5j^q 


(4.26) 

(4.27) 

(4.28) 

(4.29) 


which  in  compact  form  is 


Z(t)  -  A  (t)Z(t)  +  HC-b.u 
n  3  f 


(4.30) 


where  Z(t)  -  and 


®  ■  [•  o1  • 


The  homogeneous  part  of 


(4.30)  is  uniformly  asymptotically  stable  for  u  sufficiently  rich.  Thus  Z(t) 
is  bounded  but  not  zero  since  u  is  a  bounded  general  input  signal.  The  loss 
of  robustness  Is  due  to  the  fact  that  the  adaptive  observer  is  designed  for 
strictly  proper  systems  and  therefore  the  signal  036 ^u  in  the  proper  system 
(4 .13) , (4 .14)  cannot  be  seen  by  the  algorithm. 


In  the  following  section  we  show  that  the  robustness  of  the  adaptive 


observer  can  be  established  by  using  a  first  order  low  pass  filter  at  the 
output  of  the  plant. 


4.4.  Filter  Design 

The  strong  observability  of  the  fast  parasitic  modes  is  the  cause 
for  the  loss  of  robustness  of  the  adaptive  observer.  We  will  show  that  this 
can  be  avoided  at  the  cost  of  increasing  the  order  of  the  adaptive  observer 
and  adaptive  laws  by  one.  Consider  the  following  first  order  filter  for 
filtering  the  output  of  the  unknown  plant 


w«-aw  +  dy  ,a>0 

a  S''  *3 

y  ■  w. 

•'s 

Then  (4.31),  (4.32)  augmented  with  (4.5)-(4.7)  becomes 


(4.31) 

(4.32) 


2  ■  Az  +  bu  +  Hn 
o  o 

Jin  -  A^n  +  ubjU 

y  ■  [1  0  ...  0]z  ■  w 

•^s  '  o 


(4.33) 

(4.34) 

(4.35) 


where  A’ 


d  !  d  C 
si  so 

.  b« 

.Yl 

A* 

0  !  A 
i  0 

b 

o 

,  H- 

.*12  . 

and  “  (WjX^)"^.  The 


only  restrictions  on  a^,  d^  is  that  controllability  and  observability  of  the 
augmented  system  should  be  preserved.  Thus  without  loss  of  generality  we  can 
express  (4.33),  (4.35)  in  the  "modal"  canonical  form 
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Mn  ■  A^n  +  wb^u 

■  [10  . . .  Ojz  “  w. 


(4.37) 


(4.38) 


The  order  of  (4.36)  Is  n+l  and  A  and  ere  defined  appropriately  as  In  (4.9). 

The  (n-4-l)-th  adaptive  observer  for  estimating  z  and  Identifying  a  and  b  under 

S  8 

the  assumption  that  n*0  Is  described  by  (4.15)  through  (4.20)  where  y  is  replaced 

by  y  and  the  signals  In  the  algorithm  are  of  order  n+l. 
s 

The  stability  of  the  (n+l)-th  order  adaptive  observer  In  the  presence 
of  parasltlcs  Is  equivalent  to  the  stability  of  the  following  set  of  differential 


equations 


e  -  Ke  +  d[4.^v+'?^q]  -  Hn 


"i  “  ®i"y-ys 

♦  ■  -rtj^v 


(4.39) 

(4.40) 

(4.41) 

(4.42) 


which  In  compact  form  become 


Z(t)  -  A  (t)Z(t)  +  H  n 

n  n 


where  Z(t)  and  A  (t)  are  defined  appropriately  as  In  (4.25)  and  H 


(4.43) 


’  -h' 

'  "  -  -  . 
.  0. 


following  theorem  establishes  the  robustness  of  the  (n+l)-th  order  adaptive 
observer  with  respect  to  parasltlcs  and  gives  a  bound  for  the  composite  error. 
Theorem  4.4.1;  If  u(t)  is  dominantly  rich  for  the  plant  (4.36)  to  (4.38), 
then  the  homogeneous  part  of  (4.43)  is  uniformly  asymptotically  stable  and 
the  composite  error  vector  Z(t)  Is  bounded  by 

“l  “l 

Urn  sup|Z(t)l  <  UY  IHUb-l  (4.44) 

t-*<D  ^2  ®2  ^ 


where  positive  constants  In  the  bound  |ifr(t,T)|  i  nj^exp(-n2(t-T))  for 
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1 


the  transition  matrix  ^(t,T)  of  (4.42).  The  constants  y,  a^,  and  are  as 
defined  In  Theorem  4.3.1. 

Proof:  The  proof  follows  directly  from  Lemma  2.1.1  and  Theorem  2.3.1  of 
Section  2  by  noticing  the  similarity  of  (4.43)  with  (2.13). 

4.5.  Discrete-Time  Identifier 

Consider  system  (3.1),  (3.2)  of  Section  3  but  with  an  output  of 

the  form 

y(k)  -  C^x(k)  +  C22(k)  (4.45) 

Instead  of  (3.3).  By  using  the  transformation 

n(k)  ■  Z(k)  +  Px(k)-B^u(k)  +  Bj(u(k)-u(k-l))  (4.46) 

and  assuming  that  the  dominant  part  is  In  the  observable  canonical  form  we 
obtain  the  following  representation  for  the  plant. 

X  (k+1)  -  A  x  (k)  +  (aC,  +  uH)n(k)  +b  u(k)  +  (aC,  +  yA. -)B-u(k-l)  (4.47) 

S  XS  Z  S  4  12  I 

n(k+l)  -  yA^  (k)  +  yA^B^[(u(k+l)-u(k))-(u(k)-u(k-l)) ]  (4.48) 

y(k)  -  [1  0  ...  0]x  (k)  (4.49) 

s 

where 


bg  -  b-yaC2A^B^  +  [B^(I-2yA^)’‘^C^  j  0]^ 

B^  -  (I-yA^)"^(B2  +  PB^),  B^  -  (I-yA^)'^^  (4.50) 

and  P  satisfies  (3.8). 
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Due  to  this  observable  canonical  form  we  can  write  (4.47)  to 


(4.49)  as 


n  n+1 

y(k)  -  ^E^aj^y(k-l)  +  ^E^b^u(k-l)  +  y^(k) 


(4.51) 


where 


y^(k)  ■  j^^j^8in(k-l) 


(4.52) 


and  Is  the  1th  row  of  aC2'*'MH. 

The  reduced-order  system  obtained  by  setting  ii*0  In  (4.51),  (4.52) 


Is 


_  n  tt+l_ 

y(k)  -  Z  i  y(k-l)  +  E  b.u(k-l). 
1*1  ^  1*1  ^ 


(4.53) 


We  point  out  that  In  Section  2  we  assume  y(k)  ■  C^x(k)  +  yC2Z(k) ,  therefore  the 

parasltlcs  were  weakly  observable.  Here  we  allow  y  (k)  ■  (k)  +  02*  (k) . 

This  means  that  as  y-^-O  a  throughput  Is  Introduced  In  the  output  which  causes 

an  extra  delayed  Input  term  b^|  j^u(k-n-l)  to  appear  In  both  the  dominant  part 

of  the  plant  and  the  reduced  order  system.  The  adaptive  Identifiers 

described  by  (3.14)  to  (3.21)  and  (3.44)  to  (3.50)  cannot  see  the  term 

b  .,u(k-n-l).  Thus  If  we  apply  the  nth  order  identifier  described  by  (3.14) 
n+l 

to  (3.21)  or  (3.44)  to  (3.50)  to  the  plant  (4.51),  (4.52)  we  will  obtain  an 
error  system  equation  of  the  form  (3.26)  or  (3.53)  but  with  an  additional 
disturbance  Input  due  to  the  unmodeled  term  b^^^u(k-n-l) .  This  disturbance 
causes  an  output/parameter  error  which  is  non-zero  even  when  the  parasltlcs 
are  Infinitely  fast,  l.e.  y^O  and  therefore  the  scheme  Is  no  longer  robust. 

The  robustness  of  discrete-time  Identifiers  with  respect  to  strongly  observ¬ 
able  parasltlcs  can  be  established  by  modifying  the  identifiers  to  take  Into 

account  the  term  b  .,u(k-n-l)  Introduced  by  the  parasltlcs.  For  the  nth-order 
nrl 


parallel  Identifier  described  by  (3.14)  to  (3.21),  we  Introduce  the  following 
modification:  The  estimation  or  adjustable  model  Is  modified  to 
n  n+l-^ 

yp(k)  -  ^E^a^(k)yp(k-i)  +  ^E^b^(k)u(k-i)  -  p^(k)0^(k-l)  (4.54) 

yp(k)  -  p^(k-l)0^(k-l)  (4.55) 

where 

P^Oc)  -  [ij(k)  ...  bj^(k)  ...  b_j(k),b^j(k)], 

e^(b-l,  .  [y_,(k-l)  ...  y_(k-n)„(k-l)  ...  „(k-n)  ,u(k-„-l)  ]".  (4.56) 

The  adaptation  laws  for  updating  p  (k)  have  the  same  form  as  (3.18)  to  (3.20) 

o 

but  Pp(k)  and  6p(k)  are  replaced  with  ^^(k)  and  ^^(k)  respectively.  The 
dimension  of  the  adjustable  gain  Fp(k)  is  defined  appropriately  and  the 
stability  condition  (3.21)  remains  the  same.  A  similar  modification  can  be 
Introduced  for  the  series-parallel  identifier  (3.44)  to  (3.50).  Applying 
the  modified  nth  order  adaptive  identifiers  to  the  plant  (4.51),  (4.52)  with 
strongly  observable  parasltlcs  and  following  the  procedure  of  Section  3,  we 
can  show  that  the  scheme  Is  robust,  that  is,  the  composite  output /parameter 
error  Is  of  0(ij). 


4.6.  Discussion  and  Example 

The  loss  of  robustness  In  the  case  of  strongly  observable  parasltlcs 
Is  due  to  the  fact  that  the  parasltlcs  introduce  a  throughput  In  the  reduced- 
order  (p-0)  system.  Continuous-time  adaptive  schemes  are  designed  for 
strictly  proper  systems  and  are  no  longer  uniformly  asymptotically  stable 
when  applied  to  proper  systems.  In  the  case  of  discrete-time  Identifiers 
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the  parasitics  Introduce  an  extra  delayed  Input  in  the  ARMA  model  which 
cannot  be  seen  by  the  identifier.  This  input  acts  as  a  persistent  distur¬ 
bance  in  the  composite  error  equation  for  all  p.  By  modifying  the  adaptive 
schemes  at  the  expense  of  increasing  their  order  by  one  we  re-establish 
robustness  with  respect  to  strongly  observable  parasitics. 

We  demonstrate  the  effect  of  strongly  observable  parasitics 
on  the  continuous-time  adaptive  observer  by  digital  simulation  for  the  plant. 


(4.57) 


px  ■  -4x,  -  2u  (4.58) 

f  r 

y  -  [1  Olx  +  Xf  (^-59) 

which  is  transformed  into 


(4.60) 


pa  ■  -4a  +  0.5u  (4.61) 

y  «  y.  (4.62) 

T  T 

The  vectors  to  be  identified  are  a*  (5,  10]  and  b*  [0.75,  -3]  .  The  observa¬ 
tion  ana  identification  results  are  summarized  in  Figs.  4.1  and  4.2.  Figures 
4.1a,b,c  show  that  although  the  parasitics  are  infinitely  fast  (i.e.  p“0)  the 
identification  and  observation  error  is  significant.  When  modeling  error  is 
present  (p-0.1)  Figs.  4.2a,b,c  show  that  the  identification  and  observation 
errors  are  more  oscillatory.  In  both  cases  (p-0,  0.1)  the  identification  results 


are  erroneous  and  the  observation  error 


considerable.  The  effect  of  weakly 


Fig.  4.4  Identification  results  for  y,  *  O.I. 
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observable  parasltlcs  on  the  Identification  and  observation  results  has  been 
simulated  In  Section  2  where  it  Is  shown  that  the  identification  error  is 
almost  linear  with  y  and  the  adaptive  results  are  exact  when  there  Is  no 
modeling  error  (y-0) .  Thus  the  introduction  of  the  low-pass  filter  at  the 
output  makes  the  parasltlcs  weakly  observable  and  gives  satisfactory  Identi¬ 
fication  results  provided  the  parasltlcs  are  sufficiently  fast  as  It  was 
emphasized  In  Section  2.  The  results  obtained  In  this  section  can  be 
extended  to  all  other  continuous  and  discrete-time  Identifiers  and  adaptive 


observers. 


5.  REDUCED-ORDER  ADAPTIVE  CONTROL 


5.1.  Introduction 

Several  attempts  have  been  made  to  formulate  and  analyze  reduced- 
order  adaptive  control  schemes.  In  [IS]  local  stability  has  been  proved  for 
a  reduced-order  Indirect  adaptive  regulator.  Efforts  on  reduced-order  direct 
adaptive  control  [22-24]  have  been  restricted  to  single  first  or  second  order 
examples  rather  than  the  general  problem.  In  these  examples  It  was  shown  by 
simulations  [22], [23],  or  "linearization"'  [24]  that  unmodeled  parasltlcs  can 
lead  to  Instability.  The  analysis  In  Section  2  of  the  effects  of  high 
frequency  plant  Inputs  on  the  performance  of  Identifiers  and  adaptive  observers 
with  parasltlcs  has  determined  that  the  Inputs  should  be  restricted  to 
dominantly  rich  Inputs.  As  a  design  concept  the  dominant  richness  requires 
that  In  the  presence  of  parasltlcs  the  richness  condition  be  satisfied 
outside  the  parasitic  range.  It  excludes  wideband  Inputs  such  as  noise  and 
square  waves  as  undesirable.  The  situation  In  adaptive  control  Is  more 
difficult  because  the  plant  Input  Is  generated  by  adaptive  feedback  which 
Incorporates  the  unknown  plant  with  parasltlcs.  The  schemes  so  far  do  not 
contain  a  mechanism  to  restrict  the  frequency  content  of  the  plant  Input.  The 
lack  of  this  mechanism  has  caused  the  loss  of  robustness  reported  In  [22-24]. 

The  two  main  results  of  this  section  are,  first,  an  estimate  of  the 
region  of  attraction  for  adaptive  regulation  and  second,  a  modification  of 
the  adaptive  laws  to  guarantee  boundedness  In  the  case  of  tracking.  The 
frequency  content  and  magnitude  of  the  reference  Input  signal,  the  speed  ratio 
M  of  slow  verstis  fast  phenomena,  the  adaptive  gain  and  Initial  conditions  are 
shown  to  have  crucial  effects  on  the  stability  of  the  adaptive  control  schemes. 
These  results  are  first  analytical  conditions  for  robustness  of  direct 


adaptive  control  with  respect  to  high  frequency  dynamics.  They  are  obtained 
for  a  continuous-time  SISO  adaptive  control  scheme  [7].  The  same  methodology 
can  be  extended  to  more  complicated  continuous  and  discrete-time  adaptive 
control  problems.  The  section  Is  organized  In  two  main  subsections.  The  first 
contains  a  simple  motivating  scalar  example  which  Illustrates  the  salient 
features  of  the  general  methodology  developed  in  the  second  subsection. 


5.2.  The  Scalar  Reduced-Order  Adaptive  Control  Problem 

We  start  with  a  simple  example  of  reduced-order  adaptive  control  In 
which  the  output  of  a  second  order  plant 

y  "ay  +  2z-u,  a  >0  (5.1) 

■^P  P  P  P 

Vii  ■  “Z  +  u  (5.2) 


with  unknown  constant  parameters  a^  and  u> 
of  a  first  order  model 


Is  required  to  track  the  state  y 

m 


y  "-ay  +  r(t)  a  >0  (5.3) 

•'m  m'^m  m 

where  u  is  the  control  input  and  r*r(t)  is  a  reference  input,  a  uniformly 
bounded  function  of  time.  This  example  serves  as  a  motivation  for  and  an 
Introduction  to  the  general  methodology  to  be  developed  in  the  next  subsection 
As  in  Section  2,  the  model-plant  mismatch  is  due  to  some  "parasitic" 
time  constants  which  appear  as  multiples  of  a  singular  perturbation  parameter 
y  and  introduce  the  "parasitic"  state  In  (5.1),  (5.2)  the  parasitic  state 


is  defined  as  n  *  z-u  resulting  Into  the  following  representation 
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un  ■  -n  -  yu  '  .5) 

where  the  "dominant"  part  (5.4)  and  "parasitic"  part  (5.5)  of  the  plant  appear 
explicitly. 

If  we  apply  to  the  plant  with  parasltlcs  (5.4),  (5.5)  the  same 
adaptive  law  which  we  would  have  applied  to  the  plant  without  parasltlcs,  that 
Is  If  we  use  the  control 

u  -  -K(t)y^  +  r(t)  (5.6) 

P 

and  the  adaptive  law 

K  -  Yeyp  Y>0  (5.7) 

we  obtain 


where 


e  -  -a  e-  (K(t)-K*)(e  +  y  )  +  2n 
m  m 

un  -  -n  +  y(Y«(e'»y„)^-K(K-a  )(e  +  y  )  +  2Kn  +  Kr-f] 

m  pm 

K  ■  Ye(e  +  y  ) 

n 


P  w 


*  A 

K*  •  a  +a  . 
m  p 


(5.8) 

(5.9) 

(5.10) 

(5.11) 


The  existing  theory  of  adaptive  control  [33,39]  guarantees  stability  properties 
for  the  case  without  parasltlcs,  y"0,  when  (5.8),  (5.9),  and  (5.10)  reduce  to 


e  -  -a  i-  (K(t)-K*)(e  +  y  ) 

HI  (u 


K  -  Ye(e  +  yjj)- 


(5.12) 

(5.13) 


Lgggnf  g»2J.;  For  any  bounded  initial  conditions  e(0),  K(0)  the  solution  e(t), 
K(t)  of  (5.12) ,  (5.13)  is  uniformly  bounded  and  lim  e(t)“0,  lim  K(t)  ■  K  where 

t-*<»  t-^*  ®  • 


66 


constant  K  is  in  general  a  function  of  e(0),  K(0) .  Furthermore,  if  r(t)  is 
s 

sufficiently  rich  then  lim  K(t)  *  K*,  independent  of  e(0) ,  K(0) . 

t  * 

The  above  example  illustrates  some  of  the  robustness  questions  to  be 
answered  in  this  section.  Given  that  the  adaptive  system  without  parasitics,  in 
this  case  (5.12),  (5.13),  possesses  properties  such  as  in  Lemma  5.2.1,  how  will 
these  properties  be  altered  by  the  parasitics  that  is,  what  are  the  stability 
properties  of  (5.8)  to  (5.10)?  Which  modification  of  the  adaptive  law  would 
help  to  preserve  some  of  the  desirable  properties?  The  perturbation  parameter 
p  provides  us  with  a  means  to  answer  such  questions  in  a  semi-quantitative  way 
using  the  orders  of  magnitude  0(u'^),  noting  that  for  p  small,  the  quantity  0(p  ) 
is  small  when  v  >  0  and  large  when  v  <  0 .  The  smallness  of  p  implies  that  the 
parasitics  are  fast  and  that  by  neglecting  them,  p-0,  we  concentrate  on  the 
slow,  that  is  the  "dominant,"  part  of  the  plant. 

As  we  shall  see  a  first  property  to  be  lost  due  to  parasitics  is 
global  stability.  In  the  case  of  regulation,  that  is  when  r(t)  ■  0,  the 

boundedness  of  the  solutions  e(t),  K(t)  and  the  convergence  of  e(t)  to  zero  as 
t-*-<»  is  preserved,  but  is  not  global.  It  possesses  a  domain  of  attraction  whose 
size  we  describe  by  estimating  the  orders  of  magi.ltudes  of  the  axes  of  an 
elipsold  Viv) .  In  the  tracking  problem,  when  r(t) 0  the  adaptive  system  with 
parasitics  such  as  (5.8)  to  (5.10)  may  not  converge  to  or  may  not  even  possess 
an  equilibrium.  A  practical  goal  is  then  to  guarantee  some  boundedness 
properties.  We  show  that  a  redesign,  which  may  sacrifice  some  properties  of 
the  ideal  system  without  parasitics,  results  in  the  convergence  from  any  point 
in  P(p)  to  a  disk  8(p)  around  the  origin  in  the  e,  n-plane.  The  design  objec¬ 
tive  is  then  to  make  P(p)  as  large  as  possible  and  B(p)  as  small  as  possible. 
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Let  us  illustrate  this  discussion  by  analyzing  the  regulation  problem  and  the 
tracking  problem  for  the  example  (5.1)  to  (5.3). 

a.  Regulation ;  In  the  regulation  problem  expressions  (5.8)  to  (5.10)  become 


r(t)  -  0,  y^(t)  -  0.  e(t)  -  y„(t)  (5.14) 

m  /  p 

y»ay+u  +  2Tl  *  (5.15) 

P  P  P 

p,!]  ■  -  Tl  -  |j,u  (5.16) 

u  -  -K(t)y  (5.17) 

P 

K  -  Yyp  (5.18) 


and  the  objective  is  to  drive  y^  to  zero  despite  the  presence  of  parasitics 
while  assuring  that  all  the  signals  in  the  closed  loop  system  (5.15)  to  (5.18) 
remain  bounded.  It  is  Important  to  note  that  the  open  loop  system  (5.15), 
(5.16)  might  not  be  stabilizable  by  constant  gain  output  feedback  for  a  given 
value  of  p.  If  this  is  the  case  then  there  is  no  hope  that  the  adaptive 
controller  (5.17),  (5.18)  will  stabilize  the  equilibrium  of  (5.15),  (5.16). 

The  following  lemma  characterizes  parasitics  for  which  a  linear  output  stabi¬ 
lizing  feedback  law  exists. 

X,tnnift  5.2.2:  There  exists  a  >  0  and  a  constant  such  that  for  all  p€  (0,Pj^] 
the  system  (5.13),  (5.16)  with  the  feedback  law 

u  ■  -  K  y  (5.19) 

o  p 

is  an  asymptotically  stable  closed- loop  system.  Furthermore, 


(5.20) 


-  -  a  >  K  >  a 
p  o  p 


We  now  establish  the  stability  properties  of  the  adaptive  control  system  (5.15) 
to  (5.18)  for  u  < 

Theorem  5.  2.1:  There  exists  u*  <  and  positive  numbers  a<  1/2,  such  that 

for  y  6  (0,y*  ]  any  solution  y  (t),  n(t),  K(t)  of  (5.15)  to  (5.18)  starting 

P 

from  the  set 


I?(y)  -  ty  .  n  ,K;  |y  |  +  |k|  <  c^y““,  |n|  <  c^u 


-a-1/2. 


(5.22) 


is  bounded  and  y^  0,  n  0,  K(t)  -*•  constant  as  t  -+■  *. 

Proof ;  Let  K,  >  a_  be  a  constant  and  consider  the  function 


y^  (K  -  K.  )^  2 

V(yp,Tl,K)  -  -^  + - 2^-^ 


(5.23) 


Observe  that  for  each  y>0,  c  >0,a>0  the  equality 


V(yp.‘n,K)  -  cy* 


(5.24) 


defines  a  closed  surface  S(y,  a»  c)  In  R  .  The  derivative  of  V  along  the  solu¬ 
tion  of  (5.5)  to  (5.18)  Is 


»  -  -<*i-%)yJ-'^^+t‘(il  +  2yp)<vyJ+K«pyp+2'pyp-2Kyp-K^yp+2Kri) 


(5.25) 


A  detailed  analysis  of  (5.25)  shows  that  there  exist  constants  a  <  1/2,  c,  y^ 


such  that  V  ^  0  for  each  y€(0,y*]  and  all  y^,  n>  K  enclosed  in  S(y>  a>  c) . 


n 


Moreover  V  »  0  only  at  the  equilibrium  yp*0,  constant.  The  same 

analysis  shows  that  there  exist  positive  constants  c^,  c^  such  that  the  set 


P(y)  ■  {yp»n»K:  |y  |  +  (K|  <  c^p  |n|  <  C2M 


-a-1/2. 


(5.26) 


is  enclosed  by  the  surface  S(p,  a>  c)  and  any  solution  of  (5.15)  to  (5.18)  start' 
ing  from  p(y)  remains  inside  5(p,  a»  c).  Furthermore  inside  S(y.  a»  c)  V  is  a 
non-increasing  function  of  time  which  is  bounded  from  below  and  hence  converges 
to  a  finite  value  V  .  Since  v  is  bounded,  v  is  uniformly  continuous  for  all  y  , 

ca  p 

n,  K  enclosed  in  S(p*  a>  c)  and  therefore  lim  V>0  1.e.  y  .».0,n-4>0  and  K 

tn-  »  P 

constant  as  t 

Remark  5.2.1;  It  can  also  be  shown  that  increasing  adaptive  gain  y  for  a  fixed 
y  reduces  the  size  of  the  domain  P(y)  and  the  stability  properties  of  Theorem 
5.2.1  can  no  longer  be  guaranteed  if  Y>0(l/y). 

Remark  5.2.2;  As  y-^0,  domain  P(y)  becomes  the  whole  space  R  ,  that  is  the 
adaptive  regulation  problem  (5.15)  to  (5.18)  is  well  posed  with  respect  to 
parasltlcs. 

Remark  5.2.3;  Theorem  5.2.1  is  more  than  a  local  result  because  it  shows  that 

given  any  bounded  initial  condition  y  (0),  n(0),  K(0) ,  there  always  exists  y* 

P 

such  that  for  each  ye  (0,y*]  the  solution  of  (5.15)  to  (5.18)  is  bounded  and 
yp->0,  n-^O,  K-*-  constant  as  t -»•“>. 

Remark  5.2.4:  Since  Theorem  5.2.1  is  only  a  sufficient  condition  it  is  of 
Interest  to  examine  whether  the  stability  properties  of  Lemma  5.2.1  are  Indeed 
lost  for  initial  conditions  outside  the  set  (5.22).  From  Lemma  5.2.2  and  the 
fact  that  K(t)  is  nondecreasing,  it  can  be  seen  that  instability  occurs  if 


K(t  )  >--  a  . 
o  y  p 
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As  an  Illustration  of  the  stability  properties  established  by  Theorem 

5.2.1,  simulation  results  for  (5.15)  to  (5.18)  with  Sp *  4  and  different  values 

of  u,Y  and  Initial  conditions  are  plotted  In  Figs.  5.1  to  5.4.  In  addition  to 

yp(t)  also  V  with  “  7  is  plotted  against  time  to  show  whether  all  the  signals 

In  the  closed  loop  remain  bounded.  In  Figs.  5.1a,b  where  p«0.05,  y=‘5, 

y  (0)»1.0,  Ti(0)  =  1.0,  K(0)  ■  3,  the  objective  of  the  regulator  is  achieved  since 
P 

yp-»-0  and  V  is  bounded.  Increasing  FpCO)  from  1  to  2.4  and  keeping  all  the 

other  conditions  the  same  as  in  Fig.  5.1,  the  regulator  falls  its  objective  and 

y^-».oo  as  shown  in  Fig.  5.2.  With  the  same  initial  conditions  as  in  Fig.  5.1a,b 

but  with  p“0.07  Instead  of  y  =  0.05,  7^'*’°°  as  indicated  in  Fig.  5.3.  Figure  5.4 

shows  the  effect  of  Increasing  the  adaptive  gain  y.  With  the  same  initial 

conditions  as  in  Fig.  5.1a,b,  but  with  y =  30  instead  of  y = 5  regulation  fails 

and  y  -*■  « . 

P 

b.  Tracking ;  Returning  now  to  the  tracking  problem  we  note  that  for  a  general 
r(t)  0  the  system  (5.8)  to  (5>10)  need  not  possess  an  equilibrium.  The  best 
we  can  expect  to  achieve  in  this  case  is  to  guarantee  that  the  solutions 
starting  in  P(ii)  remain  bounded  and  converge  to  a  disk  B(u)  around  the  origin 
in  the  e,  n-plane. 

To  prove  such  a  result  we  modify  the  adaptive  law  (5.10)  as 

k  -  -  aK  +  Ye(e  +  y  )  (5.27) 

m 

where  a  is  a  positive  design  parameter.  In  view  of  (5.27)  the  equations 
describing  the  stability  properties  of  the  tracking  problem  in  the  presence  of 
paras 1 tics  are 

a^e  -  (K(t)  -  K*)(e  +  yj  +  2T1 


e  •  - 


(5.28) 


Fig.  5.2  Adaptive  regulation  results  for  p,  =  0.05,  y  (o)  =  2.4,  Tl(o)  =  1,  K(o)  “3 
and  Y  *  5.  P 

150 

100 

yp  50 

0 

-50 

Fig.  5.3  Adaptive  regulation  results  for  ■  0.07,  y  (c/  =  1,  Tl(o)  “  1,  K(o)  *3 


(5.29) 


uTl  -  -71+  p.[Ye(e  +  y  )‘  -  K(K-a  )(e  +  y  )  +  2KT1  +  Kr  -  r] 

in  p  in 


K  ■  -  ctK  +  veCe  +  y  ) 

in 


Theorem  5. 2 . 2 :  Let  the  reference  Input  r(t)  satisfy 

1  r(t)  I  <  ,  I  r(t)  |  <  Y  t  >  0 


(5.30) 


(5.31) 


where  r^,  r^  are  given  positive  constants.  Then  there  exist  positive  constants 
tj^,  y*,  o,  a  <  1/2,  Cj^  to  c^  such  that  for  y  S  (0,y*]  every  solution  of  (5.28) 
to  (5.30)  starting  at  t  •  0  from  the  set 


D(y)  -  {e.K.ll:  |e|  +  |k|  <  ,  |T1| 


(5.32) 


enters  the  residual  set 


V  (y)  ■  {e,n,K:  (e,n)  ^  8(y),  K  €  K} 
o 


where 


(5.33) 


-®(y)  -  |e|  +  hi  <  }  and  K  -  [K:  |k|  <  ]  (5.34) 


at  t  -  t^  and  remains  In  ^^^(y)  for  all  t  ^  fuJ^thermore  c^  >  o  >  y. 

Proof:  Choosing  the  function 


V(yp,71,K) 


(5.35) 


we  can  see  that  for  each  y>0,  c^>  0,ci>  0  the  equality 


V(e,71,K)  -  c  u' 


(5.36) 
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3 

defines  a  closed  surface  S(y,a,c^)  in  R  space.  The  derivative  of  V  along  the 
solution  of  (5.28)  to  (5.30)  Is 

2 

V  -  -  a^e^  -  2k(K-K*)  -  I  +  y(Tl  +  2e)[Ye(e+yJ^- K(K- ap)(e  +  y^) 

+  2101  +  Kr  -  2a  e  -  r  -  2 (K  -  K*)  (e  +  y^)  +  211]  (5.37) 

HI  in 

Hence 


+  |K(K-ap)(e  +  y^)|+2|KTl|  +  2|l0l|  +  2 1  (K  -  K*)  (e  +  y^^^)  | 

+  2|Tij  +  2a^(el  +  |K|rj^  +  r^}  (5.38) 

A  detailed  analysis  of  (5.38)  shows  that  there  exist  positive  constants  c  , 

0 

a<  1/2,  to  Cg,  and  y*  such  that  for  all  ye  (0,y*],  and  >  a  >  y  the  sets 

P(y) ,  ^^^(w)  given  by  (5.32)  and  (5.33),  respectively,  are  enclosed  by  S(y,a,c^) 

and  1) < 0  everywhere  Inside  S(y,a,c^),  except  possibly  in  P^(y).  Furthermore 

set  V  (y)  is  closed  and  bounded,  V  (y)cP(y)  and  P(y)/Pi  (y)  is  a  non-empty  set. 
o  o  o  . 

Thus  every  solution  of  (5.28),  (5.29),  (5.30)  starting  at  t  ■  0  from  ^^^(y)  will 

remain  in  V  (y) .  Also  every  solution  starting  at  t  ■  0  from  P(y)/0  (y)  will 
o  o 

enter  ^^^(y)  at  t  ■  t^^  and  remain  in  ^^^(y)  thereafter. 

Remark  5.2.5;  Constants  c^  i«0,l,...,4  depend  on  r^  and  which  characterize 

the  magnitude  and  frequency  content  of  the  reference  input  signal.  A  further 

analysis  of  (5.38)  indicates  that  for  a  given  y  an  Increase  in  r^^  or  can  no 

longer  guarantee  that  V<0  everywhere  in  P(y)/P  (y).  For  this  reason  our  formu- 

o 

lation  excludes  high  frequency  or  high  amplitude  reference  input  signals  such 
as  square  or  random  waveforms,  the  traditional  favorites  of  the  adaptive  control 


literature. 


Remark  5.2.6;  It  can  also  be  shown  that  increasing  the  adaptive  gain  y  for 


given  y,  r^  and  r^  reduces  the  size  of  the  domain  P(y) .  For  Y^0(l/y)  the 
stability  properties  of  Theorem  5.2.2  can  no  longer  be  guaranteed. 

Remark  5.2.7;  The  use  of  a  Is  found  to  be  essential  In  obtaining  sufficient 
conditions  for  boundedness  In  the  presence  of  parasltics.  However,  In  the 
absence  of  parasltics  (u>0) ,  a>0  causes  an  output  error  of  O(i^).  This  Is  a 
trade-off  between  boundedness  of  all  signals  In  the  presence  of  parasltics  and 
the  loss  of  exact  convergence  of  the  output  error  to  zero  In  the  absence  of 
parasltics.  The  size  of  a  reflects  our  ignorance  about  y.  If  an  upper  bound 
of  y  Is  known  a  can  be  set  equal  to  this  upper  bound.  For  high  frequency 
parasltics  y  Is  small  and  therefore  a  can  be  small. 

It  Is  of  Interest  to  examine  whether  for  Initial  conditions  outside 

the  set  P(y)  we  can  lose  boundedness.  Simulation  results  with  a  =*4,  a  »3, 

p  m 

and  y*5  are  summarized  In  Figs.  5.5  to  5.12.  Plots  of  the  output  error  and 
the  function  V  versus  time  are  obtained  for  different  Initial  conditions,  y,  a 
and  reference  Input  characteristics.  In  Fig.  5.5a,b  the  output  error  e  and 
function  V  are  plotted  against  time  for  y»0.01,  e(0)  ■  1,  n(0)  ■  1,  K(0)  »  3, 
a >0.06,  and  r(t)>3sln  2t.  The  output  error  decreases  and  remains  close  to 
zero  and  function  V  Is  strictly  decreasing  for  V  >0.05,  but  V  changes  sign 
In  the  region  V  <  0'.05  as  shown  In  Fig.  5.5b.  Keeping  the  same  conditions  as 
In  Fig.  5.Sa,b,  but  Increasing  y  from  0.01  to  0.05  we  can  still  achieve 
similar  results  as  shown  In  Fig.  5.6a,b.  However,  In  this  case  the  steady 
state  error  Is  larger  and  V  changes  sign  for  V<0.4.  Increasing  the  value  of 
y  from  0.05  to  0.08  the  output  error  becomes  unbounded  for  all  a > 0  as  Indicated 
In  Fig.  5.7.  The  effects  of  the  Input  characteristics  are  summarized  In 
Figs.  5.8  to  5.10.  In  Fig.  5.8,  y-0.05,  e(0)  -  1,  n(0)  -  1,  K(0)  -  3,  a-0,  and 


e  1000 


02  CM  QyS 

Time,t  (sec) 


Fig.  5.7  Tracking  results  for  ^  ■  0.08,  e(o)  ■  1,  TKo)  ■  1,  K(o)  ■  3,  o  >  0.0 
and  r(t)  ■  3sln2t. 


e  200 


OO  01 


02  03 

Time,!  (sec) 


Fig.  5.8  Tracking  results  for  y,  -  0.05,  e(o)  -  1,  Tl(o)  ■  1,  K(o)  -  3,  a  -  0  or  0.06 
and  r(t)  ■  38inI0t. 


Fig.  5.11  Tracking  results  for  ^  »  0,05,  e(o) 
and  r(t)  ■  3sin2t. 
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r(t)»3sinlOt  results  Into  an  unbounded  output  error  due  to  the  Increase  of 
the  frequency  of  r(t)  from  2  to  10.  The  same  Instability  result  has  been 
observed  for  a  *0.02,  0.06.  However,  for  a -0.08  the  output  error  became 
bounded  as  shown  In  Fig.  5.9  Indicating  the  beneficial  effects  of  a  when  para- 
sltlcs  are  present.  The  effect  of  the  amplitude  of  the  reference  Input  r(t)  Is 
shown  In  Fig.  5.10.  With  )i*0.05,  a*0,  0.06  and  the  same  Initial  conditions 
as  before  but  with  r(t)~15sln2t  the  output  error  goes  unbounded.  Figure  5.11 
shows  the  effect  of  Initial  conditions  on  boundedness.  By  Increasing  e(0) 
from  1  to  2.5  and  keeping  m*0.05,  n(0)  ■  1,  K(0)  ■  3,  and  r(t)  =  3sln2t  the 
output  error  becomes  unbounded  for  all  ct>0.  In  Fig.  5.12a,b  we  show  the 
loss  of  exact  convergence  of  the  output  error  to  zero  In  the  absence  of  para- 
sltlcs  (pO)  due  to  the  design  parameter  o. 

5.3.  Reduced  Order  Adaptive  Control  Problem 

We  now  consider  the  general  problem  of  adaptive  control  of  a  SISO  time* 
Invariant  plant  of  order  n  +  m  whdre  n  Is  the  order  of  the  dominant  part  of  the 
plant  and  m  Is  the  order  of  the  parasitica.  The  plant  Is  assumed  to  possess 
slow  and  fast  parts  and  Is  represented  la  the  explicit  singular  perturbation 
form 


^11*  ^12*  **1“ 

(5.40 

AjiX  +  A22*  +  b2U  . 

ReX(A22)  <  0 

(5.41 

C  X 

(5.42 

where  x,  z  are  n  and  m  vectors  respectively  and  u,  y  are  the  scalar  Input  and 
output  of  the  plant  respectively,  in  this  representation  It  Is  assumed  that  the 
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parasltlcs  are  only  weakly  observable  from  (5.42),  that  is,  the  dependence 
of  y  on  the  parasitic  modes  is  0(m).  When  parasltlcs  are  strongly  observable, 
y»Cj^x+c2Z»  constant  output  feedback  can  lead  to  instability  in  general  [40] 
and  this  case  is  of  no  Interest  to  us  at  this  moment.  However,  the  same 
filtering  technique  employed  in  Section  4.4  can  be  used  to  filter  strongly 
observable  parasitics  and  obtain  an  output  of  the  form  (5.42).  In  (5.41)  state 
z  is  formed  of  a  "fast  transient"  and  a  "quasi-steady  state"  defined  as  the 
solution  of  (5.41)  with  \iz  =  0.  This  motivates  the  definition  of  the  fast  para¬ 
sitic  state  as 


and  substituting  (5.43)  into  (5.40),  (5.41)  we  obtain  a  representation  of  (5.40), 
(5.41),  (5.42)  with  the  dominant  part  (-5.45)  and  the  parasitic  part  (5.46) 
appearing  explicitly 


• 

X  ■ 

V  '’o" 

(5.45) 

« 

m-'H  ■ 

A22^  +  p.(Aj^x  +  A^u  +  Ajl]  +  A^u) 

(5.46) 

y  • 

(5.47) 

The  output  y  of  the  system  (5.45)  to  (5.47)  is  required  to  track  the  output  y 

in 


of  an  n-th  order  reference  model 
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X  “Ax  +  b  r 

SI  SI  n  SI 


T 

y  ■  c  X 
SI  mm 


whose  transfer  function  W  (s) 

m 


W  (s)  -  c*(sl  -  A  ) 
m  m  m  m 


K 

m 


2s>(«) 

01 


(5.48) 

(5.49) 


(5.50) 


is  chosen  to  be  strictly  positive  real  and  r(t)  is  a  uniformly  bounded  reference 
input  signal. 

The  reduced-order  plant  obtained  by  setting  |j.  »  0  in  (5. 45)- (5. 47) 

Is  assumed  to  satisfy  the  following  conditions: 

(1)  The  triple  (A^,bg,c^)  completely  controllable  and  observable. 

(ii)  In  the  transfer  function 

Wo^®>  "  ®o(®^  •  V*\  ^ 

N(s)  is  a  sionic  Hurwitz  polynomial  of  degree  n-1  and  D(s)  is  a  sionic  polynomial 
of  degree  n.  For  ease  of  exposition  we  assume  that  *  ^si  ~ 

The  controller  structure  has  the  sasie  form  as  that  used  In  [7]  for  the 
parasitic  -free  plant  that  is  for  p.  -  0  in  (5.45)  to  (5.47).  In  this 
controller  the  plant  Input  u  and  measured  output  y  are  used  to  generate  a 
(2n-2)  dimensional  auxiliary  vector  v  as 

Vi  -  AVi  +  gu 


(5.52) 
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I  ^ 


-  -i 


- 


1 1 


I  i 


If 


K  c 


e 


i  -Ax  +  b  r 
m  n  in  n 


T 

y  ■  c  X 
m  mm 


(5.48) 

(5.49) 


whose  transfer  function  W  (s) 


W„(s)  -  c*(sl  -  A  ) 
in  in  IB 


-s 


K 


2n(«) 

m 


(5.50) 


Is  chosen  to  be  strictly  positive  real  and  r(t)  is  a  uniformly  bounded  reference 
input  signal. 

The  reduced -order  plant  obtained  by  setting  ^  >  0  in  (5. 45)- (5. 47) 

Is  assumed  to  satisfy  the  following  conditions: 

(1)  The  triple  (A^,b^,c^)  completely  controllable  and  observable. 

(ii)  In  the  transfer  function 

*  V*\  4 

N(s)  is  a  monic  Hurwitz  polynomial  of  degree  n-1  and  D(s)  is  a  monic  polynomial 
of  degree  n.  For  ease  of  exposition  we  assume  that 

The  controller  structure  has  the  same  form  as  that  used  in  [7]  for  the 
parasitic  -free  pUnt  that  is  for  ^  •  0  in  (5.45)  to  (5.47).  In  this 
controller  the  plant  Input  u  and  measured  output  y  are  used  to  generate  a 
(2n-2)  dimensional  auxiliary  vector  v  as 


ri  -  AVi  +  gu 


"i  -  '’"’'1 


(5.52) 


O  O  O  I 


gd*c 
®  o  o 


we  rewrite  (5.55)  to  0*36)  In  e  form  convenient  for  our  stability  analysis 


Y  -  A^Y+  b^((0  -  +  r)  + 


(5.59 


iiTl  “  ^22^  ^2®^***  ■*"  ^4®^“'*“^4^)  (5.60 

where  Aj^2  *  f^i2  ®  ^1  "  ®  ^  advantage  of  this  form  is 

that  for  d(t)  "  9*  in  the  parasitic  -free  case  (5.59)  becomes  a  non-minima I 


representation  of  the  reference  model 


*mc  “  Vmc  +  ^ 

The  equations  for  the  error  e  ■  Y  -  can  be  expressed  as 
i  -  A^.  +  b^(9  -9V(i  +  S„)  +  AjjTl 


(5.61 


(5.62 


^iTI  -  A22'n  +  M-[A^(e  +  x^)  +  A29‘’^(e  +  x^^^^)  +  A2r  +  A3:i+A^0^(e  +  Xj^) 


+  A^9^f  (9,e,Ti,r)  +  A^r] 


e,  ■  h  e  «  [1  0  .  .  0]e 


(5.63 


(5.64 


idiere 


(5.65) 


+Vi„)  +gr  +  ge^(;  +  x^) 


f(9,e.ll,r)  -  A(e^‘^+V2^)+g(«^  +  y^) 


+^12^ 


(5.66) 


^  X  •  X  j  c'  ^  A  V  “  ^  v^  “ 


(5.67) 


We  now  need  to  design  en  adepcive  law  for  updating  the  parameter  vector  6(t). 
For  the  parasltlc-free  case  [7]  the  adaptive  law 


4  ■  -  ■  -  rej^(e  +  5bc>»  r  -  >  0 


(5.68) 


guarantees  that  the  output  error  goes  to  zero  as  t  ••  •  and  the  signals  In  the 
close  loop  remain  bounded  for  any  uniformly  bounded  reference  Input  r(t).  As 
demonstrated  In  Section  5.2  for  the  scalar  tracking  problem,  the  best  we  can 
expect  in  the  presence  of  parasltics  is  to  guarantee  that  the  solutions 
starting  In  a  domain  V(,u)  remain  bounded  and  converge  to  a  disk  B(vi)  around  the 
origin  In  the  e,  n  plane.  To  achieve  this  we  modify  the  adaptive  law  (5.68)  as 


9  -  -  o9  -  raj^(«  +  x^) 


(5.69) 


where  or  Is  a  design  scalar  parameter.  The  resulting  adaptive  control  system 
with  parasltics  Is  described  by 


4-  A^e +  bg(0-0V(4  +  5^)  +1^7)  -  (5.70) 

« 

M.T|  -  AjjTl  +  ntAj(«  +  x^)  +  A20’^(S  +  ;J|^)  +  A^r  +  A^l]  -aA^0^(;  +  5;j^) 

-  A4(5  +  5a,c)^r(5  +  3^)«l+A^9^f(0,«,Tl,r)+A^tl  (5.71) 

-  [I  0  .  .  0]c  (5.72) 

e  .  -  ore  .  re^(S+5i^).  (5.73) 

Theorem  5.3.1;  Let  the  reference  Input  r(t)  satisfy 

|r(t)|<r^,  |r(t)|<r2  ^^>0  (5.74) 


for  some  given  positive  constants  r^^,  v^.  Then  there  exists  positive  constants 
y*f  a,  a<  1/2,  to  c^  and  t^^  such  that  for  each  ye  (0,y*l  every  solution  of 
(5.70)  to  (5.73)  starting  at  t  •  0  from  the  set 

r(^)  -  U.Tl.O;  l!a|l  +  ||0||  <  c^n“®,  IITIII  <  Cjy"®"''*}  (5.75) 

enters  the  resldiial  set 

(«,11)  €/9(h,),  0  (5.76) 

where 

^(H)  -  M:  llsll  +  m\  <  ,  JC  -  (9:  H  0  II  <  c^]  (5.77) 

at  t  -  t^  and  remains  inl?^(^)  for  all  t  >  tj^.  Furthermore  Cj  ><t  >  n. 


Corollary  5.3.1:  Assume  r(t)*0.  Then  there  exists  a  m*  such  that  for  all 
ye  (O.y*]  and  a»0  in  (5.73)  any  solution  e(t),  n(t),  0(t)  of  (5.70)  to 
(5.73)  which  starts  from  D(y)  given  by  (5.75)  is  bounded  and  lel-^0,  lnl-*-0, 
n  011  constant  as  t  -*■ «. 


Proof  of  Theorem  5.3.1s  Choose  the  function 


T«.“  1 


V(e,n,0)  -  ^  •  Q  )  r. .  .,C9-9-.).  +  |[T|  -  p“\e^PA^^A~^)-^]^P^[Tl 


(5. 


where  P  satisfies 


A^P  +  PA 
c  c 


■  "  -  eL 


(5. 


Pbc  ■  h  (5. 

T 

for  some  vector  q,  matrix  L  «  L  >0  and  e  >  O^and  P^  satisfies 

^1^22  ^L^l  "  ”  »  Qi  •  >  0  (5-* 

Equations  (5.79),  (5.80)  follow  from  tive  fact  that  h^(sl  -  A^)b^  is  strictly 
positive  real  [7]  and  (5.81)  follows  from  the  assumption  that  ReX(A22)  <  0. 
Observe  that  for  each  y,  >  0,  >  0,  a  <  1/2  the  equality 

V(e,Tl,0)  -  (5.i 

defines  a  closed  surface  S(y.,Qf,dQ)  in  +  derivative  of  V  along 


the  solution  of  (5.70)  to  (5.73)  is 
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V  -  -  “  e’^(qq'^  +  *L)«  -  a(9  -aVr'^e  +  p.[Tl 


-  p'^(«^PA^2A2j)^l’^Pl[Aj(e  +  x^)  +  +A2r+A3T1 


-  aA^9^(5  +  jSj^)  -A^(i  +  5^)^r(e  +  jc„^)ej^+A^9’^£(9,e,‘Il,r)+A^r 

-  P*^A2j^^P(Aee+bc(9-9V(5  +  5i^)+A^2Tl)l  (5.83) 


Let 


.-1. 


-  min  X(L),  Xj  ■  X(r  ),  X3  ■  nin  XCQ^ 


(5.84) 


Then 


V  <  -  i  Xj!|.l^,x,(l|e||  -  \\gl  I!e*||,^  -i  X3iwf  * 


+  ^IIH  -  P‘‘(  +o.j||.ll*  +03llS|l*ll.ll  +«4ll9||ll«ll 


“3ll®lr+  “6ll9||+“7l!9|lll’l!!+«*ll’l|l+ll*4ll'2+  “,1 


(5.85) 


where  to  ere  positive  constents  determined  from  r^^,  end  the  norms  of  the 
system  and  the  reference  model  matrices.  A  detailed  analysis  of  (5.85)  shows 
that  there  exists  positive  constants  9,  a  <  1/2,  Cj^  i>l,...,  5  and  m*  such 
that  for  all  y€  (0,y*]  and  c_>a>y,  p(y),  p  (y)  defined  by  (5.75)  to  (5.77) 

3  O 

are  enclosed  by  S(^3a3d^)  and  ^  <  0  everywhere  inside  except  possibly 

in  Pg(y>)<  The  set  l^(y.)  is  closed  and  bounded,  Pg(y.)  cP(y)  and  P(y)/P^(y)  is 
a  non-empty  set.  Every  solution  of  (5.70)  to  (5.73)  starting  at  t  ■  0  from 

t 

P^(m>)  will  remain  in  ^^(y.).  Since  in  P(^)/P^(y>)8  V  is  strictly  decreasing  any 
solution  starting  at  t  ■  0  from  P(y,)/P^(yk)  will  enter  at  t-  and 

remain  in  P  (y)  thereafter. 


Proof  of  Corollary  5.3.1;  The  proof  of  Corollary  5.3.1  follows  directly  from 
the  proof  of  Theorem  5.3.1  by  noting  that  when  r(t)  =  0,  x^®0,  and  a  =  0,  the 
disk  B(u)  reduces  to  the  origin  e»0,  n*0,  i.e.  in  (5.77)  c^^O. 

In  Theorem  5.3.1  and  Corollary  5.3.1  is  assumed  that  p*<  where  is 
defined  in  the  following  lemma. 

T 

Lemma  5.  3.1:  There  exists  a  y.  >  0  such  that  constant  output  feedback  u-  6  ui 
-  1  o 

stabilizes  (5.55)  to  (5.57)  for  all  ye(0,yj^]. 

The  proof  of  Lemma  5.3.1  is  more  complicated  than  that  of  Lemma  5.2.1 
and  can  be  found  in  [41]  where  an  explicit  expression  has  been  obtained  for  y^. 

Remark  5  .3.1:  The  dependence  of  constants  c^  i«l,...,4  on  r^,  shows  that 
for  a  given  y  a  reference  input  signal  with  high  magnitude  or  high  frequencies 

m 

can  no  longer  guarantee  that  V<0  everywhere  in  P(y)/P^(y).  Such  a  reference 
signal  Introduces  frequencies  in  the  input  control  signal  which  are  in  the 
parasitic  range.  Thus  the  control  signal  is  no  longer  dominantly  rich  and,  hence, 
it  excites  the  parasltlcs  considerably  and  leads  to  instability.  This  explains 
the  instability  phenomena  observed  by  other  authors  in  simulations  such  as  [23] 
where  a  square  wave  was  used  as  a  reference  input  signal. 

5.4.  Discussion 

In  this  section  we  have  analyzed  reduced-order  adaptive  control 
schemes  in  which  reference  models  can  match  the  dominant  part  of  the  plant,  while 
the  model-plant  mismatch  is  caused  by  the  neglected  high  frequency  parasitic 
modes.  In  the  presence  of  parasitics  the  global  stability  properties  of  the 
parasitic-free  schemes  can  be  lost.  However,  we  have  shown  that  in  the  regula¬ 
tion  problem  a  region  of  attraction  exists  for  exact  adaptive  regulation.  This 
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region  is  a  funcClon  of  the  adaptive  gains  and  the  speed  ratio  |i,  and  as 
H  o,  it  becomes  the  whole  space.  Thus  the  adaptive  regulation  problem  Is 
well  posed  with  respect  to  parasltlcs.  In  the  case  of  tracking  we  proposed 
a  more  robust  adaptive  law.  The  new  scheme  guarantees  the  existence  of  a 
region  of  attraction  from  which  all  signals  converge  to  a  residual  set  ^Ich 
contains  the  equilibrium  for  exact  tracking.  The  dependence  of  the  size  of 
this  set  on  design  parameters  indicates  that  a  trade-off  can  be  made 
sacrlflzlng  some  of  the  ideal  parasitic-free  properties,  in  order  to  achieve 
robustness  in  presence  of  parasitica.  The  crucial  effects  of  the  frequency 
range  of  parasltlcs,  the  adaptive  gains  and  the  reference  input  signal 
characteristics  on  the  stability  properties  of  adaptive  control  schemes, 
explain  the  undesirable  phenomena  observed  in  [22-24].  The  results  of  this 
section  are  obtained  for  a  continuous-time  SISO  adaptive  control  scheme  where 

the  transfer  ftmction  of  the  dominant  part  of  the  plant  has  a  relative  degree 

of  one.  The  same  methodology  can  be  extended  to  more  complicated  continuous 
and  discrete- time  adaptive  control  problems. 


6.  DECENTRALIZED  ADAPTIVE  CONTROL 


6.1.  Introduction  and  Problem  Statement 

In  classical  control  theory  the  systems  imder  study  are  assumed 
to  have  only  one  controller  who  determines  the  control  actions  based  on  the 
available  Information  of  the  syst^.  Systems  of  this  type  are  called 
centralized  systems. 

When  control  theory  Is  applied  to  solve  problems  of  large  systems, 
e.g.  electric  power  systems,  socioeconomic  systems,  computer  communication 
networks,  etc.,  an  Important  feature  called  decentralization  often  arises. 
Such  systems  have  several  local  control  stations;  at  each  station,  the 
controller  obseirves  only  local  system  outputs  and  controls  only  local  Inputs. 
All  the  controllers  are  Involved,  however.  In  controlling  the  overall  large 
system. 

In  decentralized  control  the  problem  of  stability  of  the  overall 
system  becomes  very  Important.  It  received  the  attention  of  many  authors  In 
the  last  few  years  [A2-4A],  but  their  studies  were  focused  mainly  on  systems 
with  known  parameters.  The  operating  environment  of  most  large-scale 
systems  Is  rather  poorly  known  and  their  parameters  cannot  be  calculated  with 
sufficient  accuracy  to  be  used  In  on-line  controllers.  Consequently,  DACRs 
(Decentralized  Adaptive  Controllers)  are  essential  for  large-scale  systems 
with  unknown  parameters. 

In  this  section  an  approach  for  the  decentralized  regulation  and 
tracking  of  a  class  of  large-scale  linear  dynamical  systems  with  unknown 
parameters  Is  developed.  Sufflclet,  condif*  s  for  decentralized  adaptive 
regulation  In  a  form  of  algebraic  cr.>.c::rla  are  established  which  can 
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guarantee  adaptive  regulation  under  certain  structural  perturbations.  In 
the  case  of  tracking,  we  propose  a  decentralized  adaptive  scheme  which 
guarantees  that  all  the  signals  of  the  closed-loop  system  are  bounded  and 
the  state  error  converges  to  a  residual  set  which  contains  the  equilibrium 
for  exact  tracking. 

The  formal  problem  statement  Is  as  follows:  Consider  a  linear 
time-invariant  system  S,  which  Is  described  as  an  Interconnection  of  N  sub¬ 
systems  S^,S2,...,S^  and  Is  represented  by 


"1 

-  A^x^  +  b^Uj^  +  f^(x) 

(6.1) 

N 

fi(x) 

-  E  A.  .X, 
j-1  J 

(6.2) 

^1 

T 

(6.3) 

1  1 

which  for  the  subsystem  :  x^€R  Is  the  state  vector,  u^€R  Is  the  control 
1  “l 

variable,  y^€  R  Is  the  output,  and  f^(x)€R  Is  the  Interaction  vector 
from  the  other  subsystems  to  S^.  The  parameters  A^,  A^^,  b^,  and  h^  are 
unknown  constant  matrices  and  all  the  triples  (A^,b^,h^)  are  completely 
controllable  and  completely  observable. 

In  this  representation  the  composite  system  S  can  be  described  as 

X  ■  Ax  +  Bu  +  f  (6.4) 

y  -  Cx  (6.5) 

T  T  T  T  T 

where  x>[x^,X2, . . .  ,x^]  Is  the  composite  state  vector,  y»  [yj^,y2»  •  •  •  »yjj] 

T 

Is  the  composite  output  vector,  u**  [u^,U2, . . . ,0^]  Is  the  composite  control 


vector,  f-Hx  represents  the  Interconnection  pattern  of  the  overall  system 
S,  and 


“  0 

^2 

^13  '  •  •  “^IN 

A„, 

0 

21 

• 

2N 

• 

• 

Al 

• 

0 

is  the  interconnection  matrix.  Furthermore,  A=diag(A^),  B>diag(b^), 
C*diag(hJ),  Ae  and  C€  R**^”  where  n*  E  n . .  The  problem  is 

to  design  LACRs  (Local  Adaptive  Controllers)  such  that  the  states  of  the 
composite  system  (6.4),  (6.5)  are  regulated  to  zero  or  track  the  states  of 
a  given  reference  model.  Each  LACR  is  associated  with  only  one  subsystem 
and  it  uses  information  only  from  that  subsystem.  That  is  the  LACRs  are 
not  allowed  to  communicate  with  each  other. 


6.2.  Adaptive  Control  Using  the  State  of  the  Plant 

To  complement  the  approach  used  in  Section  6.3  for  the  DAC 
(Decentralized  Adaptive  Control)  of  a  large-scale  system  composed  of  SISO 
subsystems  we  shall  consider  here  the  relatively  simple  case  of  a  system 
whose  state  variables  can  be  measured  and  vector  b^  is  known  for  each 
subsystem.  The  reference  model  whose  states  are  to  be  tracked  by  the 
states  of  subsystem  is 

X  -  A  X  +  b  r ,  i-l,2,...,N  (6 

m^  m^^  m^^  l  1 

"i 

where  x  €  R  ,  r .  is  a  piecewise  continuous  uniformly  bounded  reference 
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Input  and  A  Is  a  stable  matrix.  It  Is  further  assumed  that  a  vector  K* 

IIL4  X 


exists  such  that 


A.  +  b^K*  »  A 
1  11  m. 


(6.8) 


In  view  of  (6.7)  the  reference  model  M  for  the  composite  system  S  Is 


X  *  A  X  +  Br 
m  mm 


(6.9) 


T  T  T  T 

where  x  *  [x  ,x  ,...,x  ],  r ■  [r, ,r_, . . . ,r„]  ,  and  A  »  dlag(A  ).  The  basic 

™  ™1  ^  ^  ^  ^  ™  ™1 
Idea  of  the  DAC  Is  to  determine  LAC  (Local  Adaptive  Control)  Inputs  u^  such 

that  the  error  e>x-x  between  the  composite  system  S  and  reference  model  M 

ID 

as  well  as  all  signals  In  the  closed'loop  system  remain  uniformly  bounded. 

Due  to  the  Interconnection  matrix  H  it  Is  not  possible  to  Insure  11ml  el  >0 

for  any  bounded  reference  Input  vector  r.  The  maximum  we  can  achieve  In 

this  case  Is  the  convergence  of  e  to  some  bounded  residual  set.  In  the 

case  of  regulation  (r-0,  x  »0)  the  objective  is  to  find  LAC  inputs  u,  to 

m  1 

regulate  the  states  of  S  to  zero. 

a.  Regulation;  In  the  regulation  problem  r«0,  the  LAC  Inputs  u^  are 

chosen  as  [2] 


“1  “ 


(6.10) 


and  K(t)  Is  adjusted  according  to  the  adaptive  law 


*1  - 


(6.11) 


T  T 

where  r^-r^>0,  P^  =  P^>0,  and  P^  satisfies  the  Lyapunov  equation 


(6.12) 


The  closed  loop  decoupled  subsystems  are  given  by 


-  (A^  +  b^K^)x^,  i=l,2,...,N 


and  have  the  property  that  if  is  updated  according  to  (6 . 11)  x^  -»■  0 , 
K^->- constant  as  t Using  the  LACK  (6.10),  (6.11)  the  composite  system 


(6.4)  becomes 


X  -  (A  +  BK  (t))x  +  Hx 


where  K«diag(K^),  and  P“diag(P^)  is  the  solution  of 


PA^  +  A^P  -  -Q,  Q  -  diag(Q.). 
mm  1 


The  presence  of  interconnection  H  can  change  the  stability  properties  of  the 
decoupled  subsystems,  and  it  is  necessary  to  obtain  sufficient  conditions  to 
guarantee  the  stability  of  the  overall  system  S.  This  is  given  by  the 
following  theorem. 

T 

Theorem  6.2.1:  If  the  matrix  G»Q-(PH  +  H  P)  is  positive  definite  then  the 

solution  x(t),  K.  (t)  of  (6.11),  (6.14)  is  bounded  and  llmlx(t)l"0, 

1 

11ml K(t) I  "Constant. 

Proof;  Consider  the  positive  definite  function 

V  -  x^Px  +  ?  (K.-K?)’*^r7^(K  -K*).  (6.1 

X  X  X  IX 

The  time  derivative  of  V  along  the  solution  of  (6.11),  (6.14)  is 

V  -  -x^(PA  +A^P)x  +  2x’^PB(K-K*)\  +  2x’^PHx- 2  I  (K  -Kj)''^(b^P.x. )x, 
mm  i  i  i  1  1  i 


where  K*"diag(Kp.  Using  (6.15)  and  rearranging  (6.17)  we  have 
V  -  -x^(Q-  (PH+h'^P))x  -  -x^Gx. 
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If  G  Is  positive  definite  then  V  is  negative  semldeflnlte  and  therefore 

the  solution  x,  K  of  (6.1A),  (6.11)  Is  uniformly  bounded.  This  Implies  that 
»»  • 

V  is  xinlformly  bounded  and  therefore  V  is  uniformly  continuous.  Since  V  Is 
a  non-lncreaslng  function  and  Is  bounded  from  below,  It  converges  to  a  finite 
value  V  .  For  any  bounded  Initial  conditions  x(0) ,  K(0) 


11m  /  V  dt  -  V  -  V  <  « 

^  *  at  Q 


(6.19 


and  hence  llmlVI-O,  l.e.  llmixl  «0  and  llmIKI  >  constant. 
t-*«“  t-+<»  t-»<“ 

Lemma  6.2.1;  A  sufficient  condition  for  G  to  be  positive  definite  Is 


min  \  .  (Q. )  >  2|H|  max  max  A(P.) 
4  min  1  4  1 


(6.20 


Proof:  For  G  to  be  positive  definite  It  Is  sufficient  to  show  that 


x^(Q-PH-H^P)x  >  0 


(6.21 


x^Qx  >  x^(PH+h’^P)x 


(6.22 


for  all  X.  Inequality  (6.22)  Is  satisfied  If 


(6.23 


Since 


(6.24 


(6.21)  follows  from  (6.23)  and  (6.24). 

Lemma  6.2.2:  If  there  exists  a  vector  K*  such  that 


Ai  +  b^Kf 


,  1-1,2 . N 


(6.25 


where  e  Is  an  arbitrarily  small  positive  scalar  then  the  stability  properties 


of  Theorem  6.2.1  are  guaranteed  for  any  bounded  I  HI. 


Proof ;  In  view  of  (6.25)  the  properties  of  Theorem  6.2.1  are  guaranteed  If 


G  -  -  (PH+H  P) 


(6. 


Is  positive  definite.  Following  the  procedure  of  Lemna  6.2.1  G  given  by 
(6.26)  Is  positive  definite  If- 


mln  X  .  (Q. )  >  el Hi  max  max  X(P  ). 


min'^1 


£  ULLU  X  ^ 

Since  e  can  be  arbitrarily  small  (6.27)  will  always  be  satisfied  for  any 
bounded  I  Hi  . 


(6. 


b .  Tracking ;  For  the  cracking  problem  we  have  chosen  the  LAC  Inputs  u^  as 


Uj^  ■  K^(t)x^  +  r^,  i«l,2,...,N 


(6. 


and  the  adaptive  law  for  adjusting  K^(c)  as 


(6. 


where  ®  design  positive  scalar  parameter. 


Then  the 


closed-loop  composite  system  S  becomes 


X  -  A  X  +  B(K-K*)^x  +  Br  +  Hx. 
m 


(6. 


Subtracting  (6.9)  from  (6.30)  we  obtain  the  state  error  equation 


e  -  A  e  +  B(K-K*)^x  +  H(e  +  x  ) 
m  m 


(6. 


which  has  a  persistent  Input  Hx  due  to  Che  interconnections.  This  Input 
acts  as  a  disturbance  In  the  error  equation  and  therefore  the  solution  of 
(6.31)  may  not  converge  to  or  may  not  even  possess  an  equilibrium.  The 
following  theorem  establishes  sufficient  conditions  for  boundedness  and 
convergence  of  the  solution  of  (6.29),  (6.31)  to  a  residual  set. 


c 
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Theorem  6.2.2;  If  the  matrix  G*Q-(PH+H  P)  is  positive  definite  then  the 
solution  e(t),  K^(t)  1-1, 2,..., N  of  (6.29),  (6.31)  is  bounded  and  converges 

to  the  residual  sec 


,  ,2  N  7^  *  2 

P  -  =-^  Xg  +  -  SxT  >  ‘  ‘‘o> 


(6.32) 


where 


2X^1  HI  ^  2  N  *T  -1  2 

>0  - “P'V  ♦  i£i  2^  ''l  ' 


(6.33) 


X  -  min  X(G),  X.  -  min  X(r.  ),  X  -  max  max  X(P  ).  (6.34) 

g  i  1  p  i  i 

Proof:  Consider  the  po'sitive  definite  function 


T  N.  aT-I*  X 

V  -  e^Pe  +  I  (K.-K?)  r/(K.-K*) 
1-1  1111 


(6.35) 


where  P  satisfies  (6.15)  and  is  a  positive  definite  matrix.  The  time 
derivative  of  V  along  the  solution  of  (6.29),  (6.31)  is 

V  -  -e'^[Q-(PH  +  H’^P)]e  +  2e'^PHx^  +  (6.: 

T 

Assuming  that  G-Q-(PH+fl  P)  is  positive  definite  and  using  (6.34)  we  have 


V  £  -lel^X  +  2lellPllHllx  I  +  ?  ^ 

g  m  i.i  2X^  i  i 

-  2»i\(' V -- ^2xJ->  • 


Hence 


I  |2  N  2  2IPI‘IHI‘1x  I 


2,  2«  __  .2 


N  ^4  -T  _i  2 

+  Z  -rj^lK^F/l  . 
i-1  2Xj  i  i 


(6.37) 


(6.38) 
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From  (6.15)  and  (6.33) 

le|2  N  .  2 

V  <  I  o.X  (IK.I - 1-  ^  )  +d 


(6.39) 


which  Implies  that  V<0  outside  the  set  p.  Hence  any  solution  e,  K^, 
i*l,2,...,N  starting  outside  the  set  P  will  enter  P  in  finite  time  and  remain 
in  P  thereafter.  Any  solution  starting  at  t»0  from  P  will  remain  in  P  for 
all  t  >  0 . 

Remark  6.2.1;  The  results  of  Lemma  6.2.1  for  the  positive  definiteness  of 
matrix  G  applies  here  too. 

Remark  6.2.2;  If  condition  (6.25)  of  Lemma  6.2.2  is  satisfied,  we  can  show 
that  the  stability  properties  of  Theorem  6.2.2  are  guaranteed  for  any 
bounded  I  HI  . 

Remark  6.2.3;  The  use  of  is  found  to  be  essential  in  obtaining  sufficient 
conditions  for  boundedness  in  the  presence  of  unmodeled  interactions.  In 
the  absence  of  interactions,  that  is  when  each  subsystem  is  decoupled,  the 
design  parameters  > 0  cause  an  output  error  of  0(/o^).  This  is  a  trade-off 
between  boundedness  of  all  signals  in  the  presence  of  unmodeled  Interactions 
and  the  loss  of  exact  convergence  of  the  output  error  to  zero  in  the  absence 
of  interactions. 


6.3.  Adaptive  Control  Using  the  Output  of  the  Plant 

We  now  consider  the  general  problem  of  DAC  of  a  large-scale 
system  S  which  is  described  by  (6.1)  to  (6.3).  In  this  representation  only 
local  outputs  are  available  for  measurement  and  the  output  of  each  subsystem 
is  required  to  track  the  output  y^  of  an  n^-th  order  reference  model 


“i  “l  “i  ^ 
T 

C  X 


whose  transfer  function  V  (s) 

“i 

W  (s)  -  (sI-A  )"^)  -  K 

m.  HI.  m.'  m.  m.  R  (s) 

1  1  X  1  X 

Is  chosen  to  be  strictly  positive  real  and  r^(t)  Is  a  uniformly  bounded 
reference  Input  signal.  In  the  transfer  function 


T  -L  "l'*) 

W^(.)  .  h^CI-A^)  ^ 


of  the  decoupled  subsystem  It  Is  assumed  that  N^(s)  Is  a  monlc 
HuwltA  poly.<»l.l  of  iogr..  o^-l  ond  Is  .  monlc  polynonlol  of  degro. 

n. .  For  ease  of  exposition  we  also  assume  that  K . «  R  ■  l*l,2p...gN. 

X  X  m^ 

The  controller  structure  for  subsystem  Is 

+  g  u 

^  ^  ^  ^  ^  (6. 

-  Cj(t)v<^> 

.(2)  •  (2)  ^ 

vj  -  +  g^y^ 


where  Is  an  (n^-l)x(n^-l)  stable  matrix  and  (A^,gj^)  Is  a  controllable 
pair.  The  DAC  Input  Is  given  by 


u.  ■  r.  +  0.W. 


irfiere  J  ,d^(t)  ]  is  a  (2n-l) 

dimensional  adjustable  parameter  vector.  It  can  be  shown  using  the  same 
procedure  as  in  [7]  that  a  constant  vector  6*  exists  such  that  for  6^(t)  -  6* 
the  transfer  function  of  the  decoupled  subsystem  given  by  (6.43)  with 
controller  (6.44)  to  (6.46)  matches  that  of  the  reference  model  given  by 


(6.42). 


The  closed-loop  subsystem  is  described  by  the  following  equations 


0  +  g^ 


(6.47) 


where 


-  (f^(x)  0  0]^. 


(6.48) 


T  T 

Introducing  6*,  Y^«  ^ 


^i®i 


; 


Ai+8iC£  I  g^d^  , 

0  '  A, 


(6.49) 


(6.47)  is  rewritten  in  the  convenient  compact  form 


’'i  -  ''cA  *  ^1- 


(6.50) 


For  0j^«0*  and  ■  0  (6.50)  is  a  non-mlnlmal  representation  of  the  reference 
model  M. 
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The  equations  for  the  error  e.  ■  Y.-x  ,  can  be  expressed  as 

X  X  cx 


«i  -  ^ci"i  +  ^ 


(6.52) 


e^^  -  h  e^  -  [1  0  ...  0]e^. 


(6.53) 


We  now  need  to  design  an  adaptive  law  for  adjusting  the  parameter  vector 
6^(t).  For  the  decoupled  subsystem  the  adaptive  law 


0 .  *  “F  e  u 
1  1  ol  1 


Ti  -  rj  >  0 


(6.54) 


guarantees  that  e^-»-0  and  6^**-  constant  vector  as  t-^<»  for  any  uniformly 
bounded  reference  Input  r^(t).  For  the  coupled  subsystem  the  error 
equation  (6.52)  does  not  have  an  equilibrium  because  of  the  forcing  Input  F^ 
due  to  the  Interactions.  The  best  we  can  do  In  this  case  Is  to  guarantee 
boundedness  for  all  signals  In  the  closed-loop  and  convergence  of  the  state 
errors  to  a  residual  set.  We  prove  such  a  result  for  a  modified  adaptive 


0  ■  •Q  0  •  r  ©  u) 

1  11  1  ol  1 


(6.55) 


where  Is  a  positive  scalar  design  parameter.  The  state  error  equation 
for  the  DAC  scheme  Is  described  by 


e  -  A  e  +  b  (9-0*)  u  +  F(e+x  ) 
c  c  me 


(6.56) 


T 

e  ■  h^e 
o  c 


(6.57) 


where  '**  ®n^^’ 

0  -  dlag(0^),  9*  ^  diag(92),  ^  . 


A  r  T  T  T,T  A  r  iT  .  ^  fi.\ 

"  •  f"r"2 . V  *  ®o  “  f®ol’®o2 . ®oN^  ’  K  " 


In  (6.58) 


ijtj. 


The  triple  (A  ,b  ,h  )  satisfies  the  following  equations 
c  c  c 


PA  +  A  P  -  -Q  -  L 
c  c  c  c  X  c 


(6.59) 


P  b  -  h  (6.60) 

c  c  c 

where  P^- dlaR(P^^)  ,  -  dlag (q^q^)  ,  -  dlag(e^L^)  ,  ^cl“^cl^°’  L^-L^>0, 

q^  Id  a  vector  and  a  positive  scalar.  Equations  (6.59),  (6.60)  follow 

from  the  strictly  positive  realness  of  W  .(8)»h^(sI-A  .)  ^  1"1,2,...,N 

ml  cl  cl 

and  the  Kalman-Yakubovlch  lemma.  In  view  of  (6.59),  (6.60)  the  following 
theorem  establishes  the  stability  properties  of  the  DAC  scheme  for  the  overall 
system  S. 

Theorem  6.3.1:  Let 


X  ■  max  max  X(P  .),  X  ■  min  min  X(e.L, ),  X.  ■  mln(r,^). 

C£  Cl  Ej^  11  1  1 

If 

X  >  2X  IFI 
e  c 


(6.61) 

(6.62) 


then  all  the  signals  In  the  closed  loop  system  S  are  bounded.  Furthermore 
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the  solution  e(t),  6^(t),  1~1|2 . N  of  (6.55)  to  (6.57)  converges  to  the 

residual  set 


ie,2,\ 


irT^e*!  2 


V  -  (e.e^  :  -  xm)  +  ^ 


(6.63) 


where 


d  ■  E  o .  — TT - +  "  "■  -  sup 

c  i-i  1  4X.  X 

*  (f-X^IPl) 

Proof:  Choose  the  positive  definite  function 

IT  ^T— 1  j. 

V  -  f  e^P  e  +  7  E  (e.-0*)^r/(6  -0p . 
2  C  2  11  1  11 


supSx  I  . 
..  me 


(6.64) 


(6.65) 


The  time  derivative  of  V  along  the  solution  of  (6.55)  to  (6.57)  Is 


V  -  -  f  e^(Q  +L  -(P  F  +  F\))e  +  eX^x  "  E  a  (0  -0*)  ^r" V  „ 

2  C  0  C*  C  C  1DX»  <L^3»  ^  «  X  jL 


(6.66) 


Using  (6.61)  we  have 

e  2  ^  N  *  2  N 

V  <  -  X^IFI)  -  ^  I^i'^i 


2  2  2 
X^IFI^Ix 
c  me 


(rf-xm) 


(6.67) 


For  X  >  2X  IF|  V<0  outside  the  region  V  given  by  (6.63),  (6.64).  Hence 

£  C 

the  solution  0^(t),  e(t)  of  (6.55)  to  (6.57)  which  starts  outside  V  at  t*0 

will  enter  V  In  finite  time.  Once  In  V  It  cannot  escape  and  will  remain 

Inside  V  for  all  t>0.  Since  e(t)  and  x  are  bounded,  Y  and  u  are  bounded 

-  me 

and  therefore  all  the  signals  In  the  closed  loop  are  bounded. 


108 


Corollary  6.3.1;  Assume  If 

X  >  2X  |F| 
e  c 


(6.68 


and  "  0  for  1 *•  1,2, . . . ,N,  then  all  the  signals  of  the  decentralized 
adaptive  regulator  scheme  are  bounded.  Furthermore,  IYI-»-0,  I  91  ■>  constant 
as  t->-*. 

Proof:  For  x  >0  the  DAC  scheme  described  by  (6.55)  to  (6.57)  becomes  a 
— —  me 

decentralized  adaptive  regulation  problem  with  the  objective  to  regulate 

the  state  Y  to  zero.  Setting  o .  ■  0  and  x  ■  0  In  (6.55)  to  (6.57)  we  have 

1  me 


9i  - 


Y  -  A  Y  +  b  (9-9*)  u  +  FY 


Y.  -  h^Y,  Y^  ^  . 


o  c 


(6.68: 

(6.69! 

(6.70: 


Consider  the  positive  definite  function 


+  T  (e,-er)*r,  (9 -9*) 


2  iii"'!  “1'  ‘1  '“1  ''1' 


(6.7i: 


The  time  derivative  of  V  along  the  solution  of  (6.68)  to  (6.70)  Is 


V  »  -  -J-  Y^(Q  +L  -(P  F  +  f’^^P  ))Y. 
2  c  c  c  c 


(6.72: 


Using  (6.61)  we  have 

.  [  Yl  ^ 

V  <  -  (X  -  2X  IFI).  (6.73: 

2  e  c 

For  X^>2X^||F|  V<0  and  therefore  Y,  9^^  1>1,2,...,N  are  uniformly  bounded. 

••  •  • 

Since  V  Is  uniformly  bounded,  V  Is  uniformly  continuous.  Hence  11m  V»0,  l.e. 

t-«» 


SY8-*-0  and  I  91 constant  as  t-*-«. 
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6.4.  Dlsciisslon  and  Example 

In  this  section  the  problem  of  DAC  of  a  class  of  large-scale  systems 
has  been  examined.  We  have  obtained  sufficient  conditions  for  adaptive  regu¬ 
lation  and  tracking  using  DACRs.  The  proposed  approach  can  Insure  exact  con¬ 
vergence  to  the  equilibrium  point  In  the  case  of  regulation  and  boundedness 
and  convergence  to  a  residual  set  In  the  case  of  tracking.  The  size  of  the 
residual  set  depends  on  design  parameters,  the  characteristics  of  the 
reference  model  and  reference  Input  signal,  and  the  size  of  the  Interconnec¬ 
tion  matrix.  We  demonstrate  the  effectiveness  of  the  proposed  schemes  by 
digital  simulation  of  a  DAC  scheme  for  the  second  order  system 


Xj^  ■  5Xj^  +  aj^2*2  “l  (6-74) 

*2  "  ^*2  *21*1  '^2*  (6*75) 

In  this  case  we  have  two  first  order  subsystems  and  It  Is  required  to  design 
DACRs  u^  and  u^  such  that  the  states  x^,  X2  are  regulated  to  zero  or  track 


the  corresponding  states  x  , ,  x  _  of  the 

ml  m2, 

reference  model 

X  -  <■  -4x  ,  +  r, 
ml  ml  1 

(Ml) 

(6.76) 

and 

V  ■  '’*.2  *  'z 

(Mj) 

(6.77) 

respectively.  The  Interconnection  matrix  H  Is  given  by  H 
Consider  the  decoupled  subsystems 


*1  "  ^*1  “1  ^®1^ 


(6.78) 


X.  -  3x„  +  u.. 


(SJ. 


(6.79) 


a.  Regulation;  Following  the -  procedure  of  Section  6.2,  we  choose  for 

=  “^1*1*  “  ^*1  (6.81 

and  for 

“2  “  "*^2*2’  ^2  ’ 

Since  condition  (6.25)  is  satisfied  for  any  e  >  0  we  would  expect  Xj^,X2-»-0  as 
t-»-«  for  any  bounded  IHB  . 

In  Fig.  6.1a,b  the  adaptive  regulation  of  x^^,  x^  to  zero  is  shown 
for  ^  ^2l”^*  6.1c,d  we  show  the  time  response  of  the 

controller  parameters  which  converge  to  a  constant.  By  Increasing 

the  size  of  the  Interconnections  to  states  x^^,  X2  are 

still  regulated  to  zero  as  shown  in  Fig.  6.23,b.  As  a  result  of  the  stronger 
interactions  the  controller  parameters  K2  converge  to  constants  which 
are  higher  than  those  of  Fig.  6.1c,d. 

b .  Tracking ;  For  tracking  we  use  for 

Uj^  =  “  ^®l*l”'^l*^l  (6.83 


and  for  S, 


U2  »  -K2X2  +  r2. 


^2  *  ^®2’^2  ~  *^2^2 


(6.84 


,.4  A 

where  e.^  »  x-^-x^^ ,  ®2**2”*m2'  show  the  time  responses 


of  the  state  errors  e^^,  02  and  controller  parameters  K2  for 


a2j^”3,  rj^»2sint,  crj^=>0.01,  r2“3sin2t  and  a^^O.OZ.  Both  state  errors 


converge  to  a  bounded  residual  set  and  are  close  to  zero.  By  increasing 


the  interconnections  to  a^2*i0»  ®21 "  ^  leaving  all  the  other  variables 


the  same  as  in  Fig.  6.3  the  state  errors  still  converge  to  a  residual  set 


Ill 


and  the  controller  parameters  remain  bounded  as  shown  in  Fig.  6.4a,b,c,d. 

The  invariance  of  the  stability  properties  of  DACRs  (6.81)  to  (6.84)  to 
the  size  of  the  interactions  is  due  to  condition  (6.25)  of  Lemma  6.2.2,  which 
holds  for  the  system  (6 . 74) ,  (6 . 75) .  For  this  example  the  DACR  pushes  the 
controller  gains  in  a  direction  so  that  the  interactions  are  dominated.  In 
general,  however,  this  is  not  true  since  condition  (6.25)  cannot  be  satisfied 
for  all  systems  of  the  class  considered  in  this  section.  When  an  upper 
bound  for  the  norm  of  the  Interactions  is  known  the  DACRs  for  regulation  or 
tracking  can  be  designed  so  that  the  sufficient  conditions  for  stability  and 
botmdedness  obtained  in  this  section  are  satisfied.  This  can  be  achieved  by 
proper  selection  of  the  adaptive  gains  and  reference  model. 
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J. CONCLUSIONS 


In  this  thesis  we  analyze  the  robustness  properties  of  model 
reference  adaptive  schemes  with  respect  to  modeling  errors.  In  Sections  2 
to  5  the  modeling  error  consists  of  fast  parasltlcs  whereas  In  Section  6 
Is  due  to  neglected  Interactions  between  subsystems.  In  the  case  of  Identi¬ 
fiers  and  adaptive  observers  we  Introduce  the  notion  of  dominantly  rich 
Inputs,  that  Is  Inputs  which  are  rich  for  the  dominant  modes,  but  do  not 
contain  high  frequencies  In  the  parasitic  range.  In  order  to  establish 
robustness  with  respect  to  parasltlcs.  Detailed  bounds  for  the  composite 
state/parameter  error  are  established  which  Indicate  possibilities  for 
reducing  the  error  by  a  proper  choice  of  the  Input  signal.  We  show  that 
the  assumption  of  weak  observability  of  parasltlcs  In  the  plant  output  Is 
crucial  for  robustness.  When  parasltlcs  are  strongly  observable  we  show  that 
Identifiers  and  adaptive  observers  are  no  longer  robust.  In  this  case  robust¬ 
ness  Is  established  by  properly  modifying  the  adaptive  algorithms  or 
filtering  the  plant  output. 

The  effects  of  unmodeled  fast  dynamics  on  the  stability  and  perfor¬ 
mance  of  adaptive  control  schemes  are  analyzed.  We  show  that  In  the  regula¬ 
tion  problem  global  stability  properties  are  no  longer  guaranteed,  but  a 
region  of  attraction  exists  for  exact  adaptive  regulation.  In  the  case  of 
tracking  we  proposed  a  more  robust  adaptive  law  which  guarantees  the 
existence  of  a  region  of  attraction  from  which  all  signals  converge  to  a 
residual  set  which  contains  the  equilibrium  for  exact  tracking. 

The  adaptive  control  of  a  class  of  large-scale  systems  Is  also 
considered.  Decentralized  adaptive  controllers  are  proposed  for  adaptive 
regulation  and  tracking.  We  obtain  sufficient  conditions  In  the  form  of 
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algebraic  criteria  which  can  gtiaranCee  stability  or  boundedness  under 
certain  structural  perturbations.  The  theoretical  results  obtained  are 
Illustrated  by  simulations  of  numerical  examples. 

There  are  several  possible  directions  for  future  research.  The 
results  obtained  In  this  thesis  for  deterministic  algorithms  can  be  extended 
to  algorithms  In  a  stochastic  environment.  A  wider  class  of  algorithms  can 
be  analyzed  with  respect  to  more  general  structural  perturbations  following 
the  same  procedure  as  In  this  thesis.  The  stability  properties  of  Indirect 
adaptive  controllers  such  as  self-tuning  regulators  with  respect  to  modeling 
errors  such  as  fast  parasltlcs  can  be  examined.  The  decentralized  adaptive 
Identification  or  control  of  a  wider  class  of  large-scale  systems  can  be 
analyzed. 
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APPENDIX  A.  ERROR  SYSTEM  OF  MINIMAL-FORM  ADAPTIVE  OBSERVERS 


Case  1;  Adaptive  observer  f33] 

Without  loss  of  generality  let  us  assume  that  the  model  of  the 
dominant  part  of  the  plant  (2.4)  is  in  the  observable  canonical  form 


I  I 

X  =  -a  I”  ~  x  +  bu+Hn 
I  0 


(A.l) 


un  =  A^n  +  uA^^b^u 


(A.  2) 


y  =  c  X  =  [1  0  ...  0]x  =  x^. 


(A.  3) 


The  algorithm  [33]  for  the  nth  order  adaptive  observer  based  on  the 
dominant  part  (A.l),  (A. 3)  without  the  parasitics  (n = 0  in  (A.l))  is  given 
by  the  equations  (A. 4)  through  (A. 11),  below.  The  observer  equation  is 


z  =  Kz+  [k-a(t)]y  +  b(t)u  +  w+  r 


(A.  4) 


y  =  c  z  =  z, 


(A.5) 


where  w  and  r  are  auxiliary  signals  formed  by  the  output  error  e^^  =  y-y  and 
the  components 


^  _ s _  _ s _ 

,  ,  ,  *1’  '*i  “  n-1^  ,  n-2  . 

s  +  d„s  +  •  •  •  +  d  s  +  d_s  +  •  •  •  +  d 

^  n  2  n 


u  (A. 


of  the  vectors  v  and  q  as  follows 


updating  the  estimated  parameters  are  given  by 


A  ^  -V 

where  <|)  =  a-a(t)  and  'l'  =  b(t)-B  are  the  parameter  errors. 

Case  2;  Adaptive  observer  [32],  [33] 

The  following  "modal"  canonical  form  is  chosen  for  the  dominant 
part  of  the  plant  (2.4) 

X  +  bu  +  Hn  (A. 12) 


•  -1 

yn  =  A^n  +  pA^  B^u  (A. 13) 

T 

y  =  c  X  =  Xj^  (A.  14) 

T 

where  h  =  [1  1  ...  1],  A  is  an  (n-l)x(n-l)  diagonal  matrix  with  arbitrary 

but  known  constant  and  negative  diagonal  elements  (i=2,...,n)  and  a,  B 
are  the  unknown  vectors  to  be  identified.  It  is  shown  in  [5]  that  any 
completely  observable  system  can  be  represented  in  this  "modal"  canonical 
form.  The  structure  of  the  adaptive  observer  based  bn  (A. 2)  with  n = 0  is 
summarized  in  the  equations  (A. 15)  through  (A. 20).  The  adaptive  observer 
equation  is 


z  =  Kz  +  (k-a(t))y  +  b(t)u  +  w  +  r  (A. 15) 

T 

y  =  c  z  =  Zj^  (A.  16) 


where  w  and  r  are  auxiliary  signals  formed  by  the  derivatives  of  the  parameter 


error  components  and  the  components 


of  the  vectors  v  and  q  as  follows 


0 

0 

w  =  - 

• 

• 

,  r  - 

• 

d»  V 

_  n 

i  q 

L  n  J 

(A. 18) 


Moreover,  K 


1  T 
-^1  1  ^ 

I 


1  A 


L  0  I 


T  -1  1 

is  stable,  the  transfer  function  c  (sI-K)  d  =  — t-t" 

S  1  A  _ 


is  strictly  positive  real,  and  d»  (1  0  ...  0)^.  Note  that  the  first  components 
of  the  vectors  v  and  q  are  respectively.  The  adaptive  laws  for 
adjusting  the  parameters  are 


-Fe^^u  =  a(t) 


4'  *  -Mej^q  =  b(t). 


(A. 19) 
(A. 20) 


The  stability  properties  of  the  adaptive  observer  in  the  presence 
of  parasltlcs  for  Cases  1  and  2  are  described  by  the  following  error  equations 


e  ■  Ke  +  +  4'u-Hn+w  +  r,  e^^^ce 

♦  ■  FCj^u 


4*  ■  -Mej^q. 


(A. 21) 
(A. 22) 
(A. 23) 


Proposition;  For  some  vector  signals  v,  q,  w  and  r  with  v  ■  G(p)Xj^,  q  ■  G(p)u 

•  •  T'  T 

w  ■  w((^,u)  and  r  ■  r(4',q)  the  system  (A. 21)  is  input  ([x^,u,n  1  )-output  (e^^) 


equivalent  with  the  system  (A. 24)  provided  (c  ,K)  is  completely  observable 


e  -  Ke  +  d[<fr^v  +  4'^q3  -  Hr),  ■  c'^e  »  e^^  (A. 

Proof ;  From  (A. 21),  (A. 24) 

c''*'(pI-K)  +  'J'u  +  w  +  r  -  d(«f^v  +  'i’^q)]  (A. 

where  'p'  is  the  d/dt  operator.  From  (A. 25)  we  have 

where  i  denotes  the  ith  element  of  the  corresponding  vector  and  (A. 26)  is 

satisfied  by  choosing  v,  q,  w,  and  r  as  given  in  Case  1  and  Case  2.  By 

considering  the  two  equivalent  systems  (A. 21)  to  (A. 23)  and  (A. 24),  (A. 22), 

T  T  T  T 

(A.  23)  ,  it  is  obvious  that  boundedness  of  [e  ,(j»  ,?  J  will  imply  boundedness 
T  T  T 

of  [e  ,<(i  ,'t'  ].  Thus  the  stability  of  the  adaptive  observers  in  Cases  1  and 
2  in  the  presence  of  parasitics  is  equivalent  to  the  stability  of  the 
following  differential  equations 

e  =  Ke  +  d  [  (|»^v  +  "F^q  ]  -  Hn  (A .  2 


i  -  -FEj^v 
'i'  =  -Me^q 

where  K,  d,  v,  and  q  are  defined  differently  in  Cases  1  and  2. 


(A.2J 


(A.  2' 


(A. 
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APPENDIX  B.  ERROR  SYSTEM  OF  A  NONMINIMAL  FORM  ADAPTIVE  OBSERVER 
From  (2.7)  and  (2.49) 


igl  -  GCs)  -  ^ 


(B.l) 


where 


G  (s)  -  c’^(sI-A)“H 
P 


(B.2) 


Is  the  transfer  ftmctlon  of  the  plant  when  Hn  *  0.  Let 


.  n^l  I  1  n^2  •  I  . 

b,  s  +  b_s  +  ...  +  b 

■  -k  —  —  ” 

P 


-  ^  n-1  .  n-2  ^ 

s  +  a,  s  +  a.,s  +  ...  +  o 

X  A  tl 


(B.3) 


where 


a  "  Ca^»a2»  •  *  •  and  B  *  Ibj^^b^*  •  •  •  »b^) 


(B.4) 


ere  the  unknown  parameter  vectors. 


Consider  a  polynomial  w  (s  +  Xa)  which  Is  relatively  prime  to  the 

1-2 


numerator  as  well  as  the  denominator  polynomials  of  G  (s)  and  X.  ^  X.  for  1, 

P  1  j 

n 

j-2,3,...,n.  Dividing  the  numerator  and  denominator  of  G(s)  by  ir  (s  +  X^) 

1-2  i 


and  expanding  them  Into  partial  fractions  we  have 


b,  + 


1  s  +  X, 


+  ..  + 


G(s) 


,lS2 


c^ad1 (sI-A)~^H  n(s) 
u(s) 


(B.5) 


s  -  a,  - 


Note  that 


1  s  X-  *  *  *  s  +  X 
2  n 


s  -  a,  - 


1  s  +  X-  •••  s  +  X 
2  n 


cI.d1UI-A)~^  . 11  .5  Hi  S  "21 

5  ,  B  ,  ^1=2  s  +  X’  1S2  s  +  X,  ’••• 


n-1  n-2 


j^22(s+  Xi) 


i22(s  +  ^i) 


?  ] 
1-2  s  +  X^J 


^21 - 


(B.6) 


(B.5)  can  be  written  as 


y(s)  ■  j[bj^u(s)  +  aj^y(s)  +  ^^2 


n  [b^u(s)  +  a^y(s)  +  T^Hn(s)] 


s  +  X. 


(B.?: 


where 


(B.8; 


(B.7)  is  represented  In  a  block  diagonal  form  in  Figure  B.l.  The  block 

diagram  of  Figure  B.l  contains  (3n-2)  integrators  and  is  a  nonminimal  reali^ 

ration  of  the  dominant  part  of  the  plant. 

T  ~ 

The  term  h  exp(At)x(0)  in  Figure  B.l  is  added  so  that  Figure  B.l  is 
equivalent  to  the  corresponding  figure  of  a  minimal  realization  of  (B.7) 
Including  initial  conditions.  Here  x=  [x2,X2, . . . ,x^]  and  x  is  the  state  of  the 
minimal  realization  based  on  (B.7). 

A  nonminimal  state-space  representation  can  be  obtained  from 
Figure  B.l  by  defining 


l.e. 


R 


s 


Z 

s 


T  T 

[ro>ro>'*«>r  1  »  w  ■  [w*,w_,...,w  ] 
i  j  n  s  c  i  n 

(z2,Z2,...,z^1^  and  T^  -  [T^ ,T^, . . . ,T^]^ 


Hn  + 


hV*^x(0) 


(B.9 


(B.IO 


yn  *  A^n  +  mA^  b^u 


(B.l 


where  h'*’  =  [1  1 


1]  and  A»dlag(X^).  The  structure  of  the  adaptive 


observer  for  (B.IO),  (B.ll)  in  the  absence  of  parasitics  (i.e.  n  =  0,  R  *  0) 

s 

is  given  in  [32],  [33]  and  the  basic  equations  are  reviewed  below.  The 


observer  equations  are 


y  »  ®gZg  +  bj^(t)u  +  b^(t)Wg-  Xj^(y-y) 


AZ  +  hy 
s  ^ 


AW  +  hu 
s 


where  y(0)  =  0,  Z  (0)  =  0,  W  (0)  =  0.  The  adaptive  laws  for  adjusting  the 

S  8 

unknown  parameters  are  given  by 


^  -  -re^^v 


'i  =  -Me^p 


where  ()>  »  [ (aj^(t)-aj^)  ,  (a(t)-ag)^] ,  ¥  »  [ (bj^(t)-bj^)  ,  (b(t)-bg)^]  and 
aT  X  T 

[y,Z  ]  ,  p*  [u,W  ]  .  The  stability  properties  of  the  nonmlnimal  adaptive 

8  S 


observer  (B.13)  to  (B.17)  in  the  presence  of  parasitics  is  described  by 
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ti 


ii 

I 


^1  "  "^i®i  + 

^  ■  -re^v 

'F  ■  -Me^^p 

where  (B.18)  Is  obtained  by  snbstracting  (B.13)  from  (B.ll). 


(B.18) 

(B.19) 

(B.20) 


lie 

■ 

I 


APPENDIX  C.  ERROR  SYSTEM  OF  A  PARAMETERIZED  ADAPTIVE  OBSERVER 


Equation  (2. A)  can  be  represented  as 


i  »  Fx  +  gy  +  Bu  +  Hn 


(C.l) 


where 


1  0  ...  0 


-fj  0  1 


-f  0  0  ...  0 

n 


(C.2) 


and  g  satisfies  gc^  *>  A-F, 


From  (C.l) 


X  -  i^'^x  +  +  Bu(T)]dT  +/^  )l’^^*="‘^^Hn(T)dT.  (C.3) 

ft  •  * 


The  first  convolution  integral  can  be  reduced  to 


/  £^^^“‘^^[Iy(T),  Iu(T)]dT-  p*  -  M(t)p* 


(C.4) 


where  p*^  ®  [g^jB*^].  By  taking 


D(t)  -  /  £^^^’”^^Hn(T)dT 


(C.5) 


we  have 


M(t)  -  FM(t)  +  tly.Iu] 


M(0)  -  0 


(C.6) 


X(t)  -  M(t)p*  +  EXP(Ft)X(0)  +  D(t) 


(C.7) 


D(t)  -  FD(t)  +  Hn(t) 


D(0)  »  0 


(C.8) 


-1. 


un  ■  A  n(t)  +  yA  BpU 


(C.9) 
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where  p*  contains  all  the  tinknown  paraneters  of  the  matrix  A  and  vector  B 
In  (2.7).  The  adaptive  observer  [34]  for  observing  the  state  x  and 
estimating  p*  when  D(t)*0  Is  reviewed  below.  The  observer  equation  Is 


i(t)  -  M(t)p(t)  +  EXP(Ft)x(0) 


(C.ll) 


y(t)  -  c^x. 


(C.12) 


The  adaptive  law  for  updating  the  imknown  vector  p,  obtained  by 


minimizing 


.  1  2 
^1  “  2  ®1 


(C.13) 


Is  given  by  [34] 


where 


p(t)  ■  -GM  (t)cej^ 

T  A 

G  •  G  >  0  and  e^^  •  y-y . 


(C.14) 


From  (C.7),  (C.ll),  and  (C.14)  the  state  error  e»x-x  and  parameter  error 

^  ^  A 

Ap  •  p  -p  satisfy 


e  -  -M(t)GM^(t)cc^e+ [FM(t)+(Iy,Iu)]Ap  +  exp(Ft)Fe(0)-FD(t)-HTi(t)  (C.15) 


Ap  -  -QlT(t)cc^e. 


(C.16) 


APPENDIX  D.  PARALLEL  IDENTIFIER  EltROR  SYSTEM 


This  appendix  is  devoted  to  the  derivation  of  (3.26)  to  (3.28) 
from  (3.12)  and  (3.14)  to  (3.20).  Substracting  (3.14)  from  (3.12) 
yields 


y(k)  -y  (k)  -  Z  a  [y(k-l)-y  (k-1)]  +  I  [a.-a  (k)]y  (k-i) 
P  i«l  1  P  1»1  1  i  P 

+  [b,-b, (k) ]u(k-i) +  y  (k) . 

1-1  ^  ^  “ 


(D.l) 


Given  (3.24),  (3.29),  and  (3.30),  (D.l)  can  be  rewritten  as 


e(kfl)  -  Ae(k)  +  bw(k) 


(D.2) 


where 


w(k)  -  [a^-a^(k+l)  ]yg(k-i+l)  +  Ibj^-b^(k+l)  ]u(k-i+l) }  +  y^(k+l) .  (D.3) 


Given  (3.16)  and  (3.25),  (D.3)  can  be  rewritten  as 
w(k)  -  0p(k)APp(k+l)  +  y^(k+l). 


(D.4) 


Subtracting  both  sides  of  (3.16)  from  the  desired  parameter  vector  yields 


Ap  (k+1)  -  p„(k)  - 


F„(k)0^(k)v®(k+1) 
_ P  P  P 

l+0^(k)F^(k)0„(k) 
P  P  P 


(D.5) 


From  (3.14)  and  (3.18) 


p„(k-l)0„(k-l)  -  y„(k)  - 


el (k-1) F  (k-1) 0„ (k-1) v° (k) 
_£ _ 2 _ E _ 2 _ 

l+6]J(k-l)F^(k-l)0^(k-l) 
P  P  P 


(D.6) 


Substituting  (D.6)  into  (3.19)  yields 


v!(k)  -  [l+0!(k-l)F,(k-l)0^(k-l)]{y(k)-y„(k)  +  E  c  Jy(k-i)-y^(k-i)  ] } , 

P  P  i-1  ^  P 


'h 
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t:-; 


r. 


M 


I 


Using  (D.l),  (3.24),  (3.28)  and  (3.29)  converts  (D.7)  to 

v°(k)  -  [l+9^(k-l)F  (k-l)0  (k-l)][(a+c)e  (k-1) +w(k-l)l 


Coisblnlng-  (D.4)  and  (D.5)  yields 


w(k)  - 


9*(k)Ap^(k) 
P  P 


0i(k)F^(k)0„(k) 

P  P  ^ _ 


v;(k+l)  +y  (k+1). 
l+0^(k)F  (k)0^(k)  ^ 


Substituting  (D.9)  into  an  incremented  version  of  (D.8)  yields 


v°(k+l)  -  [l+0^(k)F  (k)0  (k)]((a+c)e^(k)  +  0i(k)Ap^(k)  +  y  (k+1)] 
p  PPP  PPP“ 


03(k)F„(k)0„(k)v°(k+l) 
P  P  P  P 


or 


v°(k+l)  -  (a+c)ep(k)  +  0j(k)App(k)  +  y,j(k+l) . 


Substitute  (D.4),  (D.5),  and  (D.ll)  into  (D.2)  yielding 


Cp (k+1)  -  Aep (k) +  b0p (k) App (k) +  by^(k+l) 


b0.(k)F_(k)0  (k) 


- ® - = -  [  (a+c)e  (k)  +  0  (k)  Ap  (k)  +  y  (k+1)  ] 

1+0* (k)F  (k)0  (k)  p  p  P  u 


or 


e  (fc+1) 
P 


0*(k)F„(k)0„(k)b(a+c) 
A--B - -E - B - 


l+0;(k)F^(k)0^(k) 

PPP 


ep(k)+b 


0*(k)F„(k)9„(k) 
T  _ _ 2 _ P  P 

•  M 


l+9Vk)F^(k)9^(k) 
P  P  P  . 


+  b 


0*(k)F^(k)0^(k) 

T- _ PPP 

*  <H 


l+0Vk)F„(k)0„(k) 

PPP 


y„(k+i) . 


(D.8) 

(D.9) 

(D.IO) 

(D.ll) 

(D.12) 

(k)APp(k) 

(D.13) 


similarly  substitute  (D.ll)  into  (D.5),  then 


Ap„(k+1) 

P 


App(k)  - 


F„(k)0„(k) 

P  P _ 

l+0][(k)F^(k)9^(k) 
P  P  P 


[(a+c)ep(k)  +  0p(k)APp(k)  +  y^(k+l)  1 . 

C 


After  regrouping  (D.14),  (D.13)  and  (D.14)  can  be  compactly  written  as  in 
(3.26)  to  (3.28). 
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APPENDIX  E.  SERIES-PARALLEL  IDENTIFIER  ERROR  SYSTEM 

This  appendix  derives  (3.53)  to  (3.55)  from  (3.12)  and  (3.44)  to 
(3.50).  Subtracting  (3.44)  from  (3.12)  yields 

y(k)-y3(k)  -  J^(a^-a^(k)]y(k-l)  +  Jj^[b^-b^(k)  lu(k-l)  +  y^(k)  (E.l) 

or  using  C3.46),  (3.47),  (3.52), 


y(lc+l)  -  y.(k+l)  -  0^(k)Ap,(k+l)  +  y..(W-l) 

S  9  11 


(E.2) 


where 


^Pg^k)  •  [a^  ...  a^  b^  ...  b^]  '“p^(k). 


(E.3) 


From  (3.45),  (3.49),  (3.50),  and  (E.3) 


Ap  (k+1)  •  Ap  (k)  - 

9  9 


F^(k)9(k)[e^(k)Ap^(k)  ■t»y^(fc»l)] 
l+0’^(k)F^(k)0(k) 


(E.4) 


Regrouping  (E.4)  and  substituting  It  Into  (E.2)  verifies  the  parametrlzatlon 
of  (3.53)  to  (3.55).' 
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